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Abstract 

We present a measurement of the W boson mass using data collected 
with the CDF detector during the 1994-95 collider run at the Fermilab 
Tevatron. A fit to the transverse mass spectrum of a sample of 30,115 
W ^ eu events recorded in an integrated luminosity of 84 pb~^ gives a 
mass Mw = 80.473 ± 0.065(stat.) ± 0.092(syst.) GeV/c^ A fit to the trans- 
verse mass spectrum of a sample of 14,740 W fiv events from 80 pb~^ 
gives a mass Mw = 80.465 ± 0.100(stat.) ± 0.103(syst.) GeV/c^ The dom- 
inant contributions to the systematic uncertainties are the uncertainties in 
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the electron energy scale and the muon momentum scale, 0.075 GeV/c^ and 
0.085 GeV/c^, respectively. The combined value for the electron and muon 
channel is M]y = 80.470 it 0.089 GeV/c^. When combined with previously 
pubhshed CDF measurements, we obtain Mw = 80.433 ± 0.079 GeV/c^. 

I. INTRODUCTION 

This paper describes a measurement of the W mass using W boson decays observed in 
antiproton-proton {pp) collisions produced at the Fermilab Tevatron with a center-of-mass 
energy of 1800 GeV. The results are from an analysis of the decays of the W into a muon 
and neutrino in a data sample of integrated luminosity of 80 pb^^, and the decays of the W 
into an electron and neutrino in a data sample of 84 pb~^, collected by the Collider Detector 
at Fermilab (CDF) from 1994 to 1995. This time period is referred to as Run IB whereas 
the period from 1992 and 1993 with about 20 pb^^ of integrated luminosity is referred to as 
Run lA. 

The relations among the masses and couplings of gauge bosons allow incisive tests of 
the Standard Model of the electroweak interactions |[l| . These relations include higher-order 
radiative corrections which are sensitive to the top quark mass, Mtop, and the Higgs boson 
mass, Mniggs [H- The W boson mass provides a significant test of the Standard Model 
in the context of measurements of the properties of the Z boson, measurements of atomic 
transitions, muon decay, neutrino interactions, and searches for the Higgs boson. 

Direct measurement of the W mass originated at the antiproton-proton collider at 
CERN 0. Measurements at the Fermilab Tevatron collider by CDF |^ and D0 0] have 
greatly improved precision. At LEP II, the W boson mass has been measured from the 
W pair production cross section near threshold and by direct reconstruction of the two 
Ws 0. The average of direct measurements including the analysis in this paper is of 
80.39 ±0.06 GeV/c2 §. 

Indirect W mass determinations involve Z boson measurements at LEP and SLC |Q, 
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charged- and neutral-current neutrino interactions at Fermilab [T^ , and the top quark mass 
measurement at Fermilab |1T[. A recent survey ||^ gives a W mass of 80.381 ±0.026 GeV/c^ 
inferred from indirect measurements. 

The paper is structured as follows. A description of the detector and an overview of 
the analysis are given in Section |T[ The calibration and alignment of the central track- 



ing chamber, which provides the momentum scale, is described in Section |T|. Section |I1 
also describes muon identification and the measurement of the momentum resolution. Sec- 



tion 1^ describes electron identification, the calorimeter energy scale, and the measurement 
of the energy resolution. The effects of backgrounds are described in Section 0. Section ^ 
describes a Monte Carlo simulation of W production and decay, and QED radiative cor- 



rections. Section [VI 1| describes the measurement of the detector response to the hadrons 
recoiling against the W in the event, necessary to infer the neutrino momentum scale and 
resolution. The knowledge of the lepton and recoil responses is incorporated in the Monte 
Carlo simulation of W production and decay. Section |V111| gives a description of the fitting 
method used to extract the W mass from a comparison of the data and the simulation. 
It also presents a global summary of the measured values and the experimental uncertain- 
ties. Finally, the measured W mass is compared to previous measurements and current 
predictions. 



II. OVERVIEW 

This section begins with a discussion of how the nature of W boson production and 
decay motivates the strategy used to measure the W mass. The aspects of the detector 
and triggers critical to the measurement are then described. A brief description of the data 
samples used for the calibrations and for the mass measurement follows. A summary of 
the analysis strategy and comparison of this analysis with our last analysis concludes the 
section. 
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A. Nature of W Events 



The dominant mechanism for production of W bosons in antiproton-proton coUisions is 
antiquark-quark annihilation. The W is produced with momentum relative to the center-of- 
mass of the antiproton-proton collision in the transverse {x, y) and longitudinal [z] directions 
(see Figure |I|). The transverse component of the momentum is balanced by the transverse 
momentum of hadrons produced in association with the W , referred to as the "recoil", as 
illustrated in Figure 0. 
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FIG. 1. One quarter of the CDF detector. The detector is symmetric about the interaction 
point. CDF uses a cylindrical coordinate system with the z (longitudinal) axis along the proton 
beam axis; r is the transverse coordinate, and ^ is the azimuthal angle. Pseudorapidity (ry) is 
defined as r/ = — ln(tan(^/2)), where 9 is the polar angle relative to the proton-beam direction. 
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FIG. 2. Kinematics of W boson production and decay for the events used in this analysis, as 
viewed in the plane transverse to the antiproton-proton beams. The recoil energy vector u is the 
sum of the transverse energy vectors EJj, of the particles recoiling against the W. Although energy 
is a scalar quantity, "transverse energy" commonly denotes the transverse component of the vector 
whose magnitude is the energy of the particle and direction is parallel to the momentum of the 
particle. 

The W boson decays used in this analysis are the two-body leptonic decays producing 
an electron or muon and a neutrino. Since the apparatus neither detects the neutrino nor 
measures the z-component of the recoil momentum, much of which is carried in fragments of 
the initial proton and antiproton at small angles relative to the beams, there is insufficient 
information to reconstruct the invariant mass of the W on an event-by-event basis. This 
analysis uses the transverse mass of each W event, which is analogous to the invariant mass 
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except that only the components transverse to the beamhne are used. Specifically, 

{M^f = (E^ + - (E^^ + E^)2, (1) 

where Mjf is the transverse mass of the W , is the transverse energy (see Figure 0) of the 
electron or the transverse momentum of the muon, and El^ is the transverse energy of the 
neutrino. The boldface denotes two-component vector quantities. The transverse energy of 
the neutrino is inferred from apparent energy imbalance in the calorimeters, 

?5t = E^ = -(E^ + u), (2) 

where u denotes the transverse energy vector of the recoil (see Figure ^ measured by the 
calorimeters. 



B. Detector and Triggers 

This section briefiy describes those aspects of the CDF detector and triggers pertinent 
to the W mass measurement. A more detailed detector description can be found in Refer- 
T^ ; recent detector upgrades are described in Reference and references therein. 



ence 



The CDF detector is an azimuthally and forward-backward symmetric magnetic detec- 
tor designed to study pp collisions at the Tevatron. The magnetic spectrometer consists 
of tracking devices inside a 3-m diameter, 5-m long superconducting solenoidal magnet 
which operates at 1.4 T. The calorimeter is divided into a central region (30° < 6 < 150°) 
outside the solenoidal magnet, end-plugs (10° < 9 < 30°, 150° < 9 < 170°), which form the 
pole pieces for the solenoidal magnet, and forward and backward regions (2° < < 10°, 
170° < 9 < 178°). Muon chambers are placed outside (at larger radius) of the hadronic 
calorimeters in the central region and behind added shielding. An elevation view of one 
quarter of the CDF detector is shown in Figure |1|. 
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1. Tracking Detectors 



A four-layer silicon microstrip vertex detector (SVX') |jT6[| is used in this analysis to 
provide a precision measurement of the location of the beam axis (luminous region). The 
SVX' is located directly outside the 1.9-cm radius beryllium beampipe. The four layers of 
the SVX' are at radii of 2.9, 4.3, 5.7, and 7.9 cm from the beamline. Outside the SVX' is a 
set of vertex time projection chambers (VTX) |T^, which provides r-z tracking information 
out to a radius of 22 cm for |?7| < 3.25. The VTX is used in this analysis for finding 
the z position of the antiproton-proton interaction (the event vertex). The event vertex is 
necessary for event selection, lepton track reconstruction, and the calculation of Et- 

Both the SVX' and VTX are mounted inside the central tracking chamber (CTC) ||18| , 
a 3.2-m long drift chamber that extends in radius from 31.0 cm to 132.5 cm. The CTC has 
84 sampling wire layers, organized in 5 axial and 4 stereo "super-layers" . Axial super-layers 
have 12 radially separated layers of sense wires, parallel to the z axis, that measure the r-0 
position of a track. Stereo super-layers have 6 sense wire layers, with a ~2.5° stereo angle, 
that measure a combination of r-0 and z information. The stereo angle direction alternates 
at each stereo super-layer. Axial and stereo data are combined to form a 3-dimensional 
track. Details of the calibration and alignment of the CTC are given in Section |T|. 

Track reconstruction uses r-0 information from the beam axis and the CTC axial layers, 
and z information from the VTX z vertex and the CTC stereo layers. In this analysis, the 
electron or muon momentum is measured from the curvature, azimuthal angle, and polar 
angle of the track as the particle traverses the magnetic field. 

2. Calorimeters 

The electromagnetic and hadronic calorimeters subtend 2tt in azimuth and from —4.2 
to 4.2 in pseudorapidity {r]). The calorimeters are constructed with a projective tower ge- 
ometry, with towers subtending approximately 0.1 in pseudorapidity by 15° in (central) 
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or 5° in (plug and forward). Each tower consists of an electromagnetic calorimeter fol- 
lowed by a hadronic calorimeter at larger radius. The energies of central electrons used 
in the mass measurement are measured from the electromagnetic shower produced in the 
central electromagnetic calorimeter (CEM) |T9|. The central calorimeter is constructed as 
24 "wedges" in for each half of the detector (—1.1 < f] < and < t] < 1.1). Each 
wedge has 10 electromagnetic towers, which use lead as the absorber and scintillator as the 
active medium, for a total of 480 CEM towers.0 A proportional chamber (CES) measures 
the electron shower position in the and z directions at a depth of ~ 6 radiation lengths in 



the CEM [T^. A fiducial region of uniform electromagnetic response is defined by avoiding 
the edges of the wedges. For the purposes of triggering and data sample selection, the CEM 
calibrations are derived from testbeam data taken during 1984-85; the tower gains were set 
in March 1994 using Cesium-137 gamma-ray sources. Details of the further calibration of 



the CEM are given in Section |V. 

The calorimeters measure the energy fiow of particles produced in association with the W. 
Outside the CEM is a similarly segmented hadronic calorimeter (CHA) |2^. Electromagnetic 
and hadronic calorimeters which use multi-wire proportional chambers as the active sampling 
medium extend this coverage to \ri\ = 4.2 [^. In this analysis, however, the recoil energy 
is calculated only in the region of full azimuthal symmetry, \r]\ < 3.6. Understanding the 
response of these devices to the recoil from bosons is difficult from first principles as it 
depends on details of the fiow and energy distributions of the recoil hadrons. The energy 
response to recoil energy is parameterized primarily using Z e^e~ and Z — > fi^fi^ events. 



Details of the calibration of the calorimeters to recoil energy are given in Section VII. 



* There are actually only 478 physical CEM towers; the locations of two towers are used for the 
cryogenic penetration for the magnet. 
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3. Muon Detectors 



Four-layer drift chambers, embedded in the wedge directly outside (in radius) of the 
CHA, form the central muon detection system (CMU) [Q. The CMU covers the region 
\ri\ < 0.6. Outside of these systems there is an additional absorber of 0.6 m of steel followed 
by a system of four-layer drift chambers (CMP). Approximately 84% of the solid angle for 
\r]\ < 0.6 is covered by CMU, 63% by CMP, and 53% by both. Additional four-layer muon 
chambers (CMX) with partial (70 %) azimuthal coverage subtend 0.6 < |?7| < 1. Muons from 
W decays are required in this analysis to produce a track (stub) in the CMU or CMX that 
matches a track in the CTC. The CMP is used in this measurement only in the Level 1 and 



Level 2 triggers. Details of the muon selection and reconstruction are given in Section |IT1 



4- Trigger and Data Acquisition 

The CDF trigger is a three-level system that selects events for recording to magnetic 
tape. The crossing rate of proton and antiproton bunches in the Tevatron is 286 kHz, with 
a mean interaction rate of 1.7 interactions per crossing at a luminosity of ~ 1 x 10^^ cm~^ 
sec~^, which is typical of the data presented here. The first two levels of the trigger 



consist of dedicated electronics with data paths separate from the data acquisition system. 
The third level El], which is initiated after the event information is digitized and stored. 



uses a farm of commercial computers to reconstruct events. The triggers selecting W ^ 
and W ^ HP events are described below. 

At Level 1, electrons were selected by the presence of an electromagnetic trigger-tower 
with Et above 8 GeV (one trigger tower is two physical towers, which are longitudinally 
adjacent, adjacent in pseudorapidity) . Muons were selected by the presence of a track stub 
in the CMU or CMX, and, where there is coverage, also in the CMP. 

At Level 2, electrons from W decay could satisfy one of several triggers. Some required 
a track to be found in the r-0 plane by a fast hardware processor [p^] and matched to 
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a calorimeter cluster; the most relevant required an electromagnetic cluster with Et 
above 16 GeV and a track with px above 12 GeV/c. This was complemented by a trigger 
which required an electromagnetic cluster with Et above 16 GeV matched with energy in 
the CES I^B] and net missing transverse energy in the overall calorimeter of at least 20 GeV, 



with no track requirements. The muon Level 2 trigger required a track of at least 12 GeV/c 
that matches to a CMX stub (CMX triggers), both CMU and CMP stubs (CMUP triggers), 
or a CMU stub but no CMP stub (CMNP triggers). Due to bandwidth limitations, only 
about 43% of the CMX triggers and about 39% of the CMNP triggers were recorded. 

At Level 3, reconstruction programs included three-dimensional track reconstruction. 
The muon triggers required a track with above 18 GeV/c matched with a muon stub. 
There were three relevant electron triggers. The first required an electromagnetic cluster 
with Ej- above 18 GeV matched to a track with px above 13 GeV/c with requirements 
on track and shower maximum matching, little hadronic energy behind the cluster, and 
transverse profile in z in both the towers and the CES. Because such requirements may 
create subtle biases, the second trigger required only a cluster above 22 GeV with a track 
above 13 GeV/c as well as 22 GeV net missing transverse energy in the overall calorimeter. 
The third trigger required an isolated 25 GeV cluster with no track requirement and with 
25 GeV missing transverse energy. 

Events that pass the Level 3 triggers were sorted and recorded. The integrated luminosity 
of the data sample is ~80 pb~^ in the muon sample and ~84 pb^^ in the electron sample. 

C. Data Samples 

Nine data samples are employed in this analysis. These are described briefly below and 
in more detail in subsequent sections as they are used. A hst of the samples follows: 

• The ip — s> ^^fi^ sample. A sample of ~ 500, 000 ip ~^ candidates with 2.7 < 

M^+f^- < 4.1 GeV/c^ is used to investigate the momentum scale determination and to 
understand systematic effects associated with track reconstruction. 
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The T — > /i"*"//" sample. A sample of ~ 83, 000 T — > candidates with 8.6 < 

M^i+ii- < 11-3 GeV/c^ offers checks of the momentum scale determination that are 
statistically weaker but systematically better than those from the ip — > sample. 

The Z — > Ai^/x" sample. A sample of ~ 1,900 dimuon candidates near the Z mass 
determines the momentum scale and resolution, and is used to model the response of 
the calorimeters to the recoil particles against the Z and W boson, and to derive the 
Z and W pt distributions in the W ^ iiv analysis. 

The W ^ jjiv sample. A sample of ~ 14, 700 W ^ jiv candidates is used to measure 
the W mass. 

The inclusive electron sample. A sample of ~750,000 central electron candidates 
with Et > S GeV is used to calibrate the relative response of the central electromag- 
netic calorimeter (GEM) towers. 

The Run lA inclusive electron sample. A sample of ~210,000 central electron 

candidates with Et > 9 GeV is used to measure the magnitude and the distribution of 
the material, in radiation lengths, between the interaction point and the GTG tracking 
volume. 

The W ^ ev sample. A sample of ~30,100 W ^ eu candidates is used to align the 
GTG, to compare the GEM energy scale to the momentum scale, and to measure the 
W mass. 

The Z — > e~^e~ sample. A sample of ~1,500 dielectron candidates near the Z mass 
is used to determine the electron energy scale and resolution, to model the response 
of the calorimeters to the recoil particles against the Z and W boson, and to derive 
the Z and W pt distributions in the W ^ ev analysis. 

The minimum bias sample. A total of ~ 2, 000, 000 events triggered only on a 
coincidence of two luminosity counters is used to help understand underlying event. 
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D. Strategy of the Analysis 



The determination of the momentum and energy scale^ is crucial to the W mass mea- 
surement. Momentum is the kinematic quantity measured for muons; for electrons, the 
energy measured in the calorimeter is the quantity of choice as it has better resolution and 
is much less sensitive than the momentum to the effects of bremsstrahlung The spec- 
trometer measures the momentum (p) of muons and electrons, and the calorimeter measures 
the energy {E) of electrons. This configuration allows in situ calibrations of both the mo- 
mentum and energy scales directly from the collider data. The final alignment of the CTC 
wires is done with high momentum electrons, exploiting the charge independence of the 
electromagnetic calorimeter measurement since both positives and negatives should give the 
same momentum for a given energy. The momentum scale of the magnetic spectrometer is 
then studied using the reconstructed mass of the ip and T fi'^fJ'^ resonances, ex- 

ploiting the uniformity, stability, and linearity of the magnetic spectrometer. Similar studies 
for the calorimeter are done using the average calorimeter response to electrons (both e+ 
and e~) of a given momentum. The momenta of lepton tracks from W decays reconstructed 
with the final CTC calibration typically change from the initial values used for data sample 
selection by less than 10%; their mean changes by less than 0.1%. The final CEM calibration 
differs from the initial source/testbeam calibration in early runs on average by less than 2%, 
with a gradual decline of ~5% during the data-taking period. Fits to the reconstructed 
Z yU+yU^ and Z e+e~ masses, along with linearity studies, provide the final momen- 
tum and energy scales. The mass distributions are also used to determine the momentum 
and energy resolutions. 

The detector response to the recoil u is calibrated primarily using Z —>■ fi'^fi~ and 
Z e+e~ decays in the muon and electron analyses, respectively. These are input to 

'^Throughout this paper, momentum measurements using the CTC are denoted as p, and calorime- 
ter energy measurements are denoted as E. 
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fast Monte Carlo programs which combine the production model and detector simulation. 

The observed transverse mass lineshape also depends on the transverse and longitudinal 
W momentum spectra. The pj^ spectrum is derived from the Z —>■ e^e^ and Z fJ'^f^~ 
data and the theoretical calculations. The pf spectrum is measured from the leptons in the 
Z decays by taking into account the lepton momentum and energy resolution. The theo- 
retical calculations are used to correct the difference between the pf and pj^ distributions. 
The observed u distributions provide consistency checks. The longitudinal spectrum is con- 
strained by restricting the choice of parton distribution functions (PDFs) to those consistent 
with data. 

To extract the W mass, the measured W transverse mass spectrum is fit to fast Monte 
Carlo spectra generated at a range of W masses. Electromagnetic radiative processes and 
backgrounds are included in the simulated lineshapes. The uncertainties associated with 
known systematic effects are estimated by varying the magnitude of these effects in the 
Monte Carlo simulation and refitting the data. 



E. Comparison with Run lA Analysis 

This analysis is similar to that of our last (Run lA) measurement 0], with datasets 
~ 4.5 times larger. The direct use of the Z events in modeling W production and recoil 
hadrons against the W @,|12| is replaced with a more sophisticated parameterization p8 



In this analysis our efforts to set a momentum scale using the and T dimuon masses and 
then to transfer that to an energy scale using E/p for W electrons did not produce a self- 
consistent picture, particularly the reconstructed mass of the Z with electron pairs. Instead 
we choose to normalize the electron energy and muon momentum scales to the Z mass, 
in order to minimize the systematic effects, at the cost of a modest increase in the overall 
scale uncertainty due to the limited Z statistics. A discussion of this problem is given in 
Appendix A. The instantaneous luminosity of this dataset is a factor of ~2 larger, resulting 
in higher probability of having additional interactions within the same beam crossing. Also, 
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we have included muon triggers from a wider range of polar angle. 



III. MUON MEASUREMENT 

In the muon channel, the W transverse mass depends primarily on the muon momentum 
measurement in the central tracking chamber (CTC). This section begins with a description 
of the reconstruction of charged-particle trajectories and describes the CTC calibration and 
alignment. It then describes the selection criteria to identify muons and the criteria to 
select the W ^ fiu and Z candidates. The momentum scale is set by adjusting 

the measured mass from Z decays to the world-average value of the Z mass pOj . 



The muon momentum resolution is extracted from the width of the Z —>■ fi^fi~ peak in the 
same dataset. The muon momentum scale is checked by comparing the T and ip masses 
with the world-average values. Since the average muon momentum is higher in Z decays 
than W decays, a correction would be necessary for the W mass determination if there were 
a momentum nonlinearity. Studies of the Z, T, and mass measurements indicate that the 
size of the nonlinearity is negligible. 



A. Track Reconstruction 

1. Helical Fit 

The momentum of a charged particle is determined from its trajectory in the CTC. The 
CTC is operated in a nearly (to within ~1%) uniform axial magnetic field. In a uniform 
field, charged particles follow a helical trajectory. This helix is parametrized by: curvature, 
C (inverse diameter of the circle in r-0); impact parameter, Dq (distance of closest approach 
to r = 0); 00 (azimuthal direction at the point of closest approach to r = 0); zq (the z 
position at the point of closest approach to r = 0); and cot 6, where 6 is the polar angle with 
respect to the proton direction. The helix parameters are determined taking into account 
the nonuniformities of the magnetic field using the magnetic field map. The magnetic field 
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was measured by NMR probes at two reference points on the endplates of the CTC during 
the data-taking period as shown in Figure and corrections are made on the magnetic field 
run-by-run to convert curvatures to momenta. 
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FIG. 3. Variation of the average magnetic field as a function of run number. The left side of 
the plot corresponds to January 1994 and the right side of the plot to July 1995. 

The momentum resolution is improved by a factor of ~2 by constraining tracks to origi- 
nate from the interaction point ("beam-constraint"). The z location of the interaction point 
is determined using the VTX for each event with a precision of 1 mm. The distribution of 
these interaction points has an RMS spread of 25—30 cm, depending on accelerator condi- 
tions. The r-0 location of the beam axis is measured with the SVX', as a function of 2;, to 
a precision of 10 fim. The beam axis is tilted with respect to the CTC axis by a slope that 
is typically about 400 microns per meter. 
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2. Material Effects on Helix Parameters 



The material between the interaction region and the CTC tracking volume leads to the 
helix parameters measured in the CTC that are different than those at the interaction point. 
For example, in traversing 7% of a radiation length, muons lose about 5 MeV on average 
due to dE/dx energy loss, which is significant for low tracks. Because of its small mass, 
electrons passing through the material have a large amount of (external) bremsstrahlung 
which changes both the curvature and impact parameter of the electrons. The beam con- 
straint fit accounts for the dE/dx, and restores some of the energy loss due to the external 
bremsstrahlung. In order to make accurate corrections for the dE/dx, and properly simulate 
biases from external bremsstrahlung, the magnitude and distribution of the material need 
to be understood. 

The material distribution is measured using a Run lA sample of 210,000 photon conver- 
sions, where the conversion rate is proportional to the traversed depth in radiation lengths.^ 
Conversion candidates are selected from the 9 GeV inclusive electron sample. An electron 
associated with an oppositely-charged partner track close in 6 and distance at the point 
of conversion (the point at which the two helices are parallel in azimuth) is identified as 
a 7 ^ e+e" candidate. To optimize the resolution on the measured conversion location, 
a two-constraint fit is applied to the helix parameters of the two tracks: the separation is 
constrained to vanish, and the angle from the beam spot to the conversion point is con- 
strained to match the of the photon momentum vector. These constraints give an average 
observed resolution of 0.41 cm on the conversion radius, to be compared with an expected 
resolution of 0.35 cm. The radial distributions for conversions and backgrounds up to the 
innermost superlayer in the CTC are shown in Figure |^. The prominent peak at 28 cm is due 
to the inner support structure of the CTC. Other structures such as the silicon layers of the 



•f-The Run lA and Run IB detectors are identical except for the SVX. This difference, estimated 
to be less than 0.1% of a radiation length, is neghgible compared to the total radiation length. 
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SVX and the VTX walls can be clearly resolved. This resolution is important since we need 
to fix the proportionality constant between conversions and radiation lengths by calibrating 
on a feature of known composition. The CTC inner support is chosen for this purpose since 
its construction is well-documented. Its thickness at normal incidence is (1.26 ± 0.06)% of a 
radiation length. The result for the integrated material thickness before the CTC volume, 
averaged over the vertex distribution and angular distribution, is (7.20 ± 0.38)% of a radi- 
ation length |. Variations in conversion-finding efficiency and electron trigger efficiency as 
a function of the conversion point are taken into account. Other choices for the "standard 
radiator" such as the wires of the innermost superlayer in the CTC, as shown in Figure |^, 
give consistent results. 



^This value is for electrons from W decay. Due to difference in the detector acceptance between 
electrons and muons, the material thickness for muons is (7.10 it 0.38)%. 
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FIG. 4. The radial (R) distributions for conversions (solid line) and background (dashed line) 
for the Run lA inclusive electron sample. R is negative when the photon momentum direction is 
opposite to the vector from the beam spot to the conversion position due to the detector resolution. 
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FIG. 5. Reconstructed photon conversion vertex density in the r — (f) plane for the innermost 
superlayer in the CTC, folded into 1/30 of the circumference (this layer has 30-fold symmetry). 
Each point represents one reconstructed vertex. 

Another check is provided by the E/p distribution Q of electrons from W decay (see 
Figure where E is the electron energy measured by the CEM and p is the electron 
momentum measured by the CTC. External bremsstrahlung photons pO| are collinear with 
the electron track at emission and typically point at the calorimeter tower struck by the 
electron track so that the calorimeter collects the full energy. Since the track momentum is 



** For convenience, the requisite factor of c is dropped in the ratio E/p. 
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reduced by the radiated energy, the E/p distribution develops a high-side tail. Final state 
radiation from electron production (internal bremsstrahlung) is about a 20 % contribution 
to this tail. We define the fraction of events in the tail, ftaiii to be the fraction of events 
in the region 1.4 < E/p < 1.8. The lower bound is far enough away from the peak to be 
insensitive to resolution effects. After a small QCD background correction, we find : 

f^^^ = 0.0488 ± 0.0014(stat.) ± 0.0004(syst.). 

The Monte Carlo simulation, including internal radiative effects, reproduces this value when 
the material equals (7.55 ± 0.37)% of a radiation length, in good agreement with the value 
from conversion photons above. 

An appropriate material distribution is applied to muon and electron tracks on a track- 
by-track basis. 
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FIG. 6. E/p distribution of electrons in the W ^ ev sample. The histogram indicates the 
simulation. 



B. CTC Calibration and Alignment 

The CTC calibration and alignment proceeds in two steps. First, the relationship be- 
tween the measured drift time and the distance to the sense wire is established. Second, the 
relative alignment of wires and layers in the CTC is performed. Small misalignments left 
after these procedures are removed with parametric corrections. 



1. Time-to- distance calibration 

Electronic pulsing, performed periodically during the data-taking period, gives relative 
time pedestals for each sense wire. Variations in drift properties for each super-layer are 
removed run-by-run. Additional corrections for nonuniformity in the drift trajectories are 
made based on data from many runs. After the calibration and alignment described in 
Section [ill B 2| , the CTC drift-distance resolution is determined to be 155 fim (outer layers) 



to 215 fim (inner layers), to be compared with ~ 120 /im expected from diffusion alone, and 



26 



~ 200 /iin expected from test-chamber results. 



2. Wire and layer alignment 
The initial individual wire positions are taken to be the nominal positions determined 



during the CTC construction |[T^. The distribution of differences between these nominal 
positions and the positions determined with an optical survey has an RMS of 25 fim. The 
84 layers of sense wires are azimuthally aligned relative to each other by requiring the ratio 
of energy to momentum E/p for electrons to be independent of charge. A physical model 
for these misalignments is a coherent twist of each endplate as a function of radius. A 
sample of about 40,000 electrons with 0.8 < E/p < 1.2 from the W ^ ev sample (see 
Figure ^ is used for the alignment. The alignment consists of rotating each entire layer on 
each end of the CTC by a different amount r x with respect to the outermost superlayer 
(superlayer 8) where the relative rotation of two endplates is expected to be the smallest 
according to the chamber construction. The stereo alignment is adjusted to account for the 
calculated endplate deflection due to wire tension. The measured deviation of each layer 
from its nominal position after this alignment is shown in Figure |^. 
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FIG. 7. The deviation (r x A(/>) of each CTC layer from its nominal position at the end plates 
{\z\ = 150 cm) in cm, versus the layer number. The solid (open) circles represent the west (east) 
CTC endplate. 

Figure § demonstrates the elimination of misalignment after the alignment (open circles). 
A small residual dependence of the J/ip mass on cot6' remains, which is removed with the 
correction. 



cote 1.0004 X cote. 



(3) 



The only significant remaining misalignments are an azimuthally(0)-modulated charge 
difference in < E/p > and a misalignment between the magnetic field direction and the 
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axial direction of the CTC. The (f) modulation is removed with the correction 

C^C- 0.00031 X sin (0o - 3.0) (4) 

where C equals to Q x l/px (GeV/c)~^, Q is the charge of the lepton, the coefficient 
corresponds to a nominal beam position displacement of 37 /im, and (p is in radians. The 
magnetic field misalignment is removed with the correction 

\C\ \C\ ■ (1 - 0.0017 • cot^ • sin(0o - 1.9)). (5) 
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FIG. 8. Measured J/ip mass relative to the overall final mass measurement as a function of 
Acot0 = cot9^+ — coiO^- . The solid triangles and open circles are before and after the Run IB 
calibration and alignment, respectively. Solid circles show the distribution with the cot^ correction 
of 1. 0004 X cote. 
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C. Muon Identification 



The W mass analysis uses muons traversing the central muon system (CMU) and the 
central muon extension system (CMX). 

The CMU covers the region \ri\ < 0.6. The CMX extends the coverage to jr^l < 1. 
There are approximately five to eight hadronic absorption lengths of material between the 
CTC and the muon chambers. Muon tracks are reconstructed using the drift chamber 
time-to-distance relationship in the transverse (0) direction, and charge division in the 
longitudinal (z) direction. Resolutions of 250 fim in the drift direction and 1.2 mm in z 



are determined from cosmic-ray studies p2|- Track segments consisting of hits in at least 
three layers are found separately in the r-0 and r-z planes. These two sets of segments are 
merged and a linear fit is performed to generate three-dimensional track segments ("stubs"). 
Figure |^ shows the effects of the bandwidth limitation of the CMX and CMNP triggers (see 



Section [II B 4|) and partial azimuthal coverage (see Section [II B 3| ) . 
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FIG. 9. The r) and distributions of muons are shown in (a) and (b) for W decays, and (c) and 
(d) for Z decays. Points (histograms) show the data (the simulation) with statistical uncertainties. 



Muons from W, Z, T, and ip decays are identified in the following manner. The muon 
track is extrapolated to the muon chambers through the electromagnetic and hadronic 
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calorimeters. The extrapolation must match to a track segment in the CMU or CMX. 
For high px muons from W ot Z decays, the r x matching is required to be within 2 cm; 
the RMS spread of the matching is 0.5 cm. For low pr muons from T and ip decays, a px 
dependent matching is required to allow for multiple scattering effects. Since the energy in 
the CEM tower(s) traversed by the muon is 0.3 GeV on average, the GEM energy is required 
to be less than 2 GeV for W and Z muons. This cut is not applied to muons from T or 
decays since T's and ■^'s are often produced with particles associated with the same initial 
partons. Since the energy in the GHA tower(s) traversed by the muon is 2 GeV on average, 
the GHA energy is required to be less than 6 GeV. In order to remove events with badly 
measured tracks, muon tracks are required to pass through all nine superlayers of the GTG, 
and to have the number of GTG stereo hits greater than or equal to 12. Muon tracks in the 
W ^ fiv and Z —>■ /x^/x" data samples must satisfy |-Do| < 0.2 cm, where Dq is the impact 
parameter in the r-0 plane of the muon track with respect to the beam spot. This reduces 
backgrounds from cosmic rays and QGD dijet events. Additional cosmic ray background 
events are removed from the W ^ fiu and Z — > samples when the hits of the muon 

track and the hits on the opposite side of the beam pipe, back-to-back in 0, can be fit as 
one continuous trojectory. 



D. Event Selection: W fiu; Z,T,ijj 

1. W ^ fiu and Z — ^ tJ''^fJ'~ event selection 

The event selection criteria for the W ^ /iiy mass measurement are intended to produce a 
sample with low background and with well-understood muon and neutrino kinematics. These 
criteria yield a sample that can be accurately modeled by simulation, and also preferentially 
choose those events with a good resolution for the transverse mass. The Z sample is used to 
calibrate the muon momentum scale and resolution, to model the energy recoiling against 
the Z and W, and to derive the Z and W transverse momentum spectra (pf and Pt)- In 
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order to minimize biases in these measurements, the Z event selection is chosen to 

be as similar as possible to the W —>■ fiu event selection. 

Both W ^ fiu and Z —>■ fj^^fi^ sample extractions begin with events that pass a Level 3 
high-pr muon trigger as discussed in Section 2. From these, a final sample is selected with 
the criteria listed in Table |I| and described in detail below. The event vertex chosen is the 
one reconstructed by the VTX closest in z to the origin of the muon track, and it is required 
to be within 60 cm in z of the origin of the detector coordinates. For the Z sample, the 
two muons are required to be associated either with the same vertex or with vertices within 
5 cm of each other. For the W sample, in order to reduce backgrounds from Z — > fi'^fi~ 
and cosmic rays, events containing any oppositely charged track with > 10 GeV/c and 
M^^track > 50 GeV/c^ are rejected. Candidate W ^ jjiu events are required to have a 
muon CTC track with pt > 25 GeV/c and a neutrino transverse energy E!^ > 25 GeV. A 
limit on recoil energy of |u| < 20 GeV reduces QCD background and improves transverse 
mass resolution. Candidate Z yU'^/i" events are required to have two muons with px > 
25 GeV/c. The two muon tracks must be oppositely charged. This requirement removes no 
events, indicating that the background in the Z sample is negligible. The transverse mass 
in the region 65 < Mt < 100 GeV/c^ and the mass in the region 80 < M < 100 GeV/c^ 
are used for extracting the W mass and the Z mass, respectively. These mass cuts apply 
only for mass fits and are absent when we otherwise refer to the W or Z sample. The final 
W sample contains 23,367 events, of which 14,740 events are in the region 65 < Mt < 
100 GeV/c^. The final Z sample contains 1,840 events which are used for modeling the 
recoil energy against the W and for deriving p^ , of which 1,697 events are in the region 
80 < M < 100 GeV/c^. 
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TABLES 



CrilxM'ioii 


W (nniis after cut 


Z (n'C'uts after eut 


Initial sample with Z vertex requirement 


60,607 


4,787 


£;CEM < 2 GeV 


56,489 


3,349 


Not a cosmic candidate 


42,296 


2,906 


Impact parameter -Do < 0-2 cm 


37,310 


2,952 


Track - muon stub match 


36,596 


2,752 


Stereo hits > 12 


34,062 


2,442 


Tracks through all CTC superlayers 


33,887 


1,991 


■PT > 25 GeV/c 


28,452 


1,966 


> 25 GeV 


24,881 


N/A 


|u| < 20 GeV 


23,367 


N/A 


pl^'' < 45 GeV/c, 70 < M'^'^ < 110 GeV/c^ 


N/A 


1,840 


Mass fit region 


14,740 


1,697 



TABLE I. Criteria used to select the W ^ jiv and Z ji^ ji samples. 
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2. T,tp ^ n'^/j, event selection 

Samples of T(1S, 2S, 3S) — > //"^/i" events and ■0(18, 2S) — > //''"//~ events are used to check 
the momentum scale determined hy Z ^ events. The sample extraction begins with 

events that pass a Level 2 and 3 dimuon trigger with muon pt > 2 GeV/c. The requirement 
on the event vertex is identical to that for the Z — > selection. Both muons are required 



to have opposite charges. 


Sample 


# of events 


T(1S) 




12,800 


T(2S) 




3,500 


T(3S) 




1,700 


J/0 




228,900 


0(2S) 




7,600 



TABLE II. The number of events in the T and samples after background subtraction. 
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Backgrounds are estimated from the dimuon invariant mass distributions in the side- 
bands (regions outside the mass peaks). The numbers of T and ip events after background 
subtraction are hsted in Table The average pt of muons in the T sample is 5.3 GeV/c, 
and that in the ip sample is 3.5 GeV/c. The distributions of muon px and the opening angle 
between the two muons in are shown in Figure |TU[ For comparison, the average pt of the 
muons and the average opening angle in the Z sample are 43 GeV/c and 165°, respectively. 
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FIG. 10. (a) Transverse momentum distributions of muons and (b) opening angle distributions 
between and fi~ in the T(1S) and J/^p samples. The histograms are normalized to unit area. 
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E. Event Selection Bias on Mw 



The W ^ nu selection requires muons at all three trigger levels. Of these, only the 
level-2 trigger has a significant dependence on the kinematics of the muon; its efficiency 
varies by ~5% with r] of the tracks. This variation, however, leads to a negligible variation 
(~2 MeV/c^) on the W mass since the Mt distribution is approximately invariant under 
pz boosts. The W mass would be more sensitive to the pt dependence of the inefficiency 
since Mt is directly related to pt- No pt dependence is seen, but the statistical limitation 
on measuring such a dependence leads to a 15 MeV/c^ uncertainty on the W mass. 

The muon identification requirements may also introduce a bias on the W mass. For 
example, if the W decays such that the muon travels close to the recoil, there is greater 
opportunity for the recoil particles to cause the muon identification to fail. These biases are 
investigated by tightening the muon identification requirements and measuring the subse- 
quent shifts in M^y. The maximum shift observed of 10 MeV/c^ is taken as a systematic 
uncertainty. 



F. Momentum Scale and Resolution 



A sample oi Z ^ /i+yU events is used to determine the momentum scale by normalizing 
the reconstructed Z A*^/^^ mass to the world- average mass |^|, and to measure the 



momentum resolution in the high-p^- region. Since the muon tracks from Z decays have 
curvatures comparable to those for the W mass determination, the systematic uncertainty 
from extrapolating the momentum scale from the Z mass to the W mass is small. The 
measurement is limited by the finite statistics in the Z peak. 

The Z fJ'^fJ'^ Monte Carlo events are generated at various values of Z mass with the Z 



width fixed to the world average p9[. The generation program includes the 7 — /i^/i events 



and QED radiative effects, Z — * /i/i7 ||3T| , p2| , but uses a QCD leading order calculation so 
that the Z is generated at = 0. The Z is then given a transverse momentum whose 
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spectrum is extracted from the Z data (see Section [V|) . The generated muons are 

reconstructed by the detector simulation where CTC wire hit patterns, measured from the 
real W —>■ eu data, are used to determine a covariance matrix of the muon track, and the 
track parameters are smeared according to this matrix. A beam constraint is then performed 
with the identical procedure as is used for the real data. The final covariance error matrix 
is scaled up by a free parameter to make the beam constraint momentum resolution agree 
with the data. The detector acceptance is modeled according to the nominal geometry. The 
simulation includes the effects of the bandwidth limitation of the CMX triggers. Figure ^ 
illustrates how well the effects of the acceptance and the bandwidth limitation are simulated. 



The mass distribution of the Z data, shown in Figure |TT], is then fit to simulated 

lineshapes, where the input Z mass and the scale parameter to the covariance matrix (or 
the momentum resolution) are allowed to vary. 



38 




0.998 0.999 1 1.001 1.002 

/ M/DG 



FIG. 11. Results of fit to Z mass and momentum resolution, (a) Invariant mass distribution. 
The points are the data, and the solid line is the Monte Carlo simulation (normalized to the data) 
with best fit. (b) Correlation between the scale factor and the momentum resolution. 

Fitting the invariant mass distribution in the region 80 < M^^ < 100 GeV/c^ with a 
fixed V z [H] yields 



Mz = 91.110 ± 0.097(stat.) ± 0.020(syst.) GeV/c^ (6) 
and momentum resolution 
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5{1/pt) = (0.091 ± 0.004(stat.)) x 10"^ (GeV/c)-^ (7) 



Equation ^ results in the momentum scale factor 



1.00085 ±0.00106 (8) 



which is applied to momenta of muons and electrons. The fit is shown in Figure |TT|. The 
two parameters, 6{1/pt) and M^^'^/M^^^, are largely uncorrelated, as shown. 



Table |T| contains a list of the systematic uncertainties on the Z mass. The largest uncer- 



tainty is from the radiative effects due to using the incomplete theoretical calculation ^1 



the calculation includes the final state radiation only and has a maximum of one radiated 
photon. The effect arising from the missing diagrams is evaluated by using the PHOTOS 
package which allows two photon emissions, and by using the calculation by U. Baur et 



al. [0 who have recently developed a complete 0{a) Monte Carlo program which incorpo- 
rates the initial state QED radiation from the quark-lines and the interference of the initial 
and final state radiation, and includes a correct treatment of the final state soft and virtual 
photonic corrections. When the PHOTOS package is used in the simulation instead, the 
change in the Z mass is less than 10 MeV/c^. The effect of the initial state radiation and 
the initial and final state interference is estimated to be 10 MeV/c^ P^. To be conservative 



these changes are added linearly and 20 MeV/c^ is thus included in the systematic uncer- 
tainty. The choice of parton distribution functions and that of the spectrum contribute 
negligible uncertainties. 

A number of checks are performed to ensure that these results are robust and unbiased. 
The masses and resolutions at low and high r] are measured to be consistent. The resolution 
is cross-checked using the E/p distribution in W ^ eu events, which is sensitive to the 



combined E and p resolution (see Section [IV F| and Figure [T^). Consistent results are found 



when much simpler techniques are used, that is, comparing the mean Mz, in the interval 
86 - 96 GeV/c^, between the data and the Monte Carlo simulation or fitting the invariant 
mass distribution with a Gaussian distribution. To address mis-measured tracks, a second 
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Gaussian term is added to smear track parameters for 8% of the Monte Carlo events. The 
change in Mz is neghgible. 



Effect 


Uncertainty on M| (MeV/c^) 


Statistics 


97 


Radiative corrections 


20 


Fitting 


negligible 


Parton distribution functions 


negligible 


spectrum 


negligible 


Detector acceptance, triggers 


negligible 


Total 


100 



TABLE III. Summary of uncertainties in measuring the Z mass. 
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G. Checks of Momentum Scale 



The momentum scale is checked using ijj and T masses, extracted by fitting the dimuon 
invariant mass distributions to simulated lineshapes which include QED radiative processes 
and backgrounds as shown in Figure 0. The muon momenta are corrected by the momentum 
scale factor shown in Eq. ^. The measured masses are summarized in Table 3.4. Table 3.5 
compares the measured masses with the world-average values. Within the momentum scale 
uncertainty, the agreement is very good. 



Resonance 


Mass (MeV/c2) 


T(1S) 


9464.3 ± 0.7(stat.) ± 1.6(syst.) ± lO.l(scale) 


T(2S) 


10028.1 ± 2.1(stat.) ± 1.6(syst.) ± 10.7(scale) 


T(3S) 


10358.9 ± 3.6(stat.) ± 1.6(syst.) ± ll.O(scale) 


J/ij 


3098.4 ± 0.1(stat.) ± l.l(syst.) ± 3.3(scale) 


^P{2S) 


3687.6 ± 0.5(stat.) ± l.l(syst.) ± 3.9(scale) 



TABLE IV. Measured masses of the T and ip resonances with the momentum scale correction. 
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Resonance 


World- Average Mass M^^*^ (MeV/c^) 


^CDF/^PDG _ 1 (%) 




9460.4 ± 0.2 


0.041 ± 0.018 ± 0.106 


T(2S) 


10023.30 ± 0.31 


0.048 ± 0.026 ± 0.106 


T(3S) 


10355.3 ± 0.5 


0.035 ± 0.038 ± 0.106 


J/i; 


3096.88 ± 0.04 


0.050 ±0.035 ±0.106 


V'(2S) 


3686.00 ± 0.09 


0.042 ± 0.033 ± 0.106 



TABLE V. Measured masses of the T and resonances with the momentum scale correction are 
compared to the world averages. The second uncertainty in the last column is the momentum scale 
uncertainty, and the first uncertainty includes the statistical and the other systematic uncertainties. 
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FIG. 12. The measured dimuon mass spectra near the (a) T masses, (b) J/V' mass, and (c) 
V'(2S) mass. The curves are the best fits of lineshapes from the Monte Carlo simulation. 



Source of Uncertainty 


Uncertainty on Mr (MeV/c^) 


Uncertainty on (MeV/c^) 


Muon energy loss 


1.5 


1.0 


Kinematics 


0.4 


0.1 


Momentum Resolution 


0.3 


0.1 


Non-Prompt Production 




0.3 
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Misalignment 


0.2 


0.1 


Background 


0.1 


0.1 


Time variation 


- 


- 


QED Radiative Effects 


0.4 


0.2 


Fitting Procedure, Window 






Total 


1.6 


1.1 



TABLE VI. Systematic uncertainties in T and ip mass measurements. 
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A list of the systematic uncertainties on the and T masses is given in Table 0. The 
entries in the table are described below. 



Muon Energy Loss: The momentum of each muon is corrected for energy loss in the 



material traversed by the muon as described in Section [III A2| . Uncertainties in the energy 
loss come from uncertainty in the total radiation length measurement and in material type. 
The measured T and ip masses vary by 0.8 MeV/c^ and 0.3 MeV/c^, respectively, when the 
average radiation length is changed by its uncertainty. Uncertainty due to material type is 
estimated to be 0.6 MeV/c^ per muon track. This leads to 1.1 MeV/c^ uncertainty in the 
T mass and 0.5 MeV/c^ uncertainty in the if) mass. There is a 0.8 MeV/c^ variation in the 
observed ip mass, which is not understood, when the mass is plotted as a function of the 
radiation length traversed. No statistically significant dependence (< 0.7 MeV/c^) on the 
total radiation length is observed in the T mass. These variations of 0.7 MeV/c^ in My and 
0.8 MeV/c^ in are taken as systematic uncertainties. Adding the uncertainties described 
above in quadrature, the total uncertainty is 1.5 MeV/c^ in Mx and 1.0 MeV/c^ in M^. 

Kinematics: Variation of the and distributions allowed by the data and cuts 
results in uncertainties of 0.4 MeV/c^ and 0.1 MeV/c^ in Mx and M^, respectively. 

Momentum Resolution: Variation of the momentum resolution allowed by the data 
results in uncertainties of 0.3 MeV/c^ and 0.1 MeV/c^ in Mx and M^, respectively. 

Non-Prompt Production: About 20% of t/^'s come from decays of B mesons, which 
decay at some distance from the primary vertex. The measured peak may be shifted by 
the application of the beam constraint. The difference in the ip mass between a fit using the 
beam constraint and a fit using a constraint that the two muons originate from the same 
vertex point is 0.3 MeV/c^. This difference is taken as an uncertainty. 

Misalignment: The CTC alignment eliminates most of the effects. The residual effects are 



measured by ip and W samples and are removed by corrections as described in Section IIIB 
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The corrections and corresponding mass shifts on My are summarized in Table [VI I| . The 
overall effects of 0.17 MeV/c^ in M-r and less than 0.1 MeV/c^ in are taken as a 
systematic uncertainty. 



Source 


Correction Formula 


AMt (MeV/c^) 


B-field direction 


|C| ^ \C\ • (1 - 0.0017 • cote • sin(</.o - 1.9)) 


+0.01 


(pQ dependence 


C - 0.00031 • sin(0o - 3.0) 


-0.24 


cote dependence 


cote 1.0004 • cote 


+0.40 


Total correction 




+0.17 



TABLE VII. Systematic uncertainties in T and tjj mass measurements. 
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Background: The backgrounds in the T and ■?/' mass peak regions are estimated by fitting 
the invariant mass distributions in the sideband regions (regions away from the peaks) 
with quadratic, hnear and exponential distributions. The backgrounds are included in the 
templates used to fit the masses. By varying the background shape, changes by less 
than 0.1 MeV/c^ and Mr changes by 0.1 MeV/c^. 

Time Variation: As shown in Figure 0, there is no indication of a time variation in the 
measured mass over the data-taking period, even though the resolution worsens due to high 
occupancy in the CTC at high instantaneous luminosity during the latter portion of the 
data-taking period. 
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FIG. 13. Variation of the measured T(1S) mass (Top) and width (Bottom) as a function of 
time. The left side of the plot corresponds to January 1994 and the right side of the plot to July 
1995. AMx is difference between the measured mass for a given time period and the mass using 
all the data. 



QED Radiative Effects: The Monte Carlo program includes final state QED radiation 
from muons. The systematic uncertainties of 0.4 MeV/c^ in Mx and 0.2 MeV/c^ in 
represent missing diagrams such as two photon emission and the interference between the 
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initial and final state radiation. 



Fitting Procedure, Window: Consistent results are found when fitting windows are 
varied or much simpler fitting techniques are used, that is, comparing the mean Mx and 
and comparing the fit results with Gaussian plus linear distributions between the data 
and the Monte Carlo simulation. 



H. Momentum Nonlinearity 

The average pt for Z decay muons is about 4.5 GeV/c higher than that for W decay 
muons. Since the momentum is calibrated with the Z mass, any nonlinearity in the mo- 
mentum measurement would translate into an incorrect momentum scale for the W mass 
measurement. The momentum nonlinearity is studied using measured masses from a wide 
range of curvatures — the CTC does not directly measure momentum, but curvature, which 
is proportional to 1/pt- The curvature ranges from 0.1 to 0.5 (GeV/c)~^ in the J/i) data, 
from 0.1 to 0.3 (GeV/c)-i in the T(1S) data, and 0.02 to 0.04 (GeV/c)-^ in the Z data. 



Figure shows the ratio of the measured mass to the world-average value as a function of 
the average curvature of two muons from these data. The ratios are fiat and all are well 
within statistical uncertainty of the ratio from the Z data. Since the curvature difference 
0.003 (GeV/c)-^ between the W and Z muons is much smaller than the range of curvature 
available in the if), T, and Z data, the nonlinearity effect in extrapolating from the Z muon 
momentum to the W muon momentum is estimated to be negligible. 
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FIG. 14. The ratio of the measured mass to the world-average value as a function of the average 
curvature or inverse momentum for the Z, T, and ip data. 

I. Summary 

The muon momentum scale is determined by normalizing the measured Z mass to the 
world-average mass. The scale in the data needs to be corrected by a factor of 1.00084 ± 
0.00106, the accuracy of which is limited by the finite statistics in the Z peak. When the 
momentum scale is varied over its uncertainty in the simulation, the measured W mass 
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changes by ±85 MeV/c^. The scale is cross-checked by and My. The momentum 
resolution, S{1/pt) = (0.091 ± 0.004) x lO^^ (GeV/c)-S is measured from the width of 
the Z —>■ 11'^ ii~ peak in the same dataset. Lepton momenta in the Monte Carlo events are 
smeared according to this resolution. When the momentum resolution is varied over its 
uncertainty in the simulation, the measured W mass changes by 20 MeV/c^. Systematic 
uncertainties due to the triggers and the muon identification requirements are estimated to 
be 15 MeV/c^ and 10 MeV/c^, respectively. 

IV. ELECTRON MEASUREMENT 

This section begins with a description of the algorithm that associates calorimeter tower 
responses with electron energy. It then describes the CEM relative cahbration procedure 
to correct for nonuniformity of the calorimeter response and time dependence. We dis- 
cuss the selection criteria to identify electrons and the criteria to select the W ^ ev and 
Z — > e^e' candidates. The electron energy scale is set by adjusting the reconstructed mass 
in Z — > e^e" decays to the world-average value of the Z mass. The electron resolution is 
measured from the width of the Z mass distribution. The electron energy scale determined 
by using the E/p distribution is discussed. A small calorimeter nonlinearity is observed, 
and a correction is applied to the electron energy for the W mass measurement. 

A. Electron Reconstruction 

The scintillation light for each tower in the CEM is viewed by two phototubes, viewing 
light collected on each azimuthal side. The geometric mean of the two phototube charges, 
multiplied by an initial calibration, gives the tower energy. For electron candidates, the 
clustering algorithm finds a CEM "seed" tower with transverse energy above 5 GeV. The 
seed tower and the two adjacent towers in pseudorapidity form a cluster. One adjacent tower 
is not included if it lies on the opposite side of the z = boundary from the seed tower. The 
total Et in the hadronic towers just behind the CEM cluster must be less than 12.5% of 
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the CEM cluster Et- The initial estimate of the electron energy is taken as the sum of the 
three (or two) CEM tower energies in the cluster. There must be at least one CTC track 
that points to the CEM cluster. The electron direction, used in the calculations of Et and 
the invariant mass, is defined by the highest pt track. The W and Z electron samples are 
further purified with additional cuts as discussed below in Section p.V C] . 

B. Uniformity Corrections 

To improve the CEM resolution, corrections are applied for known variations in response 
of the towers, dependence on shower position within the tower, and time variations over the 
course of the data-taking period. For the present measurement, the nominal uniformity cor- 
rections (testbeam) are refined using two datasets - the W electrons and the high- statistics 
inclusive electron dataset. The reference for correcting the electron energy is the track mo- 
mentum as measured by the CTC. Uniformity is achieved by adjusting the tower energy 
response (gain) until the mean E/p is flat as a function of time and (p, and agrees with the 
Monte Carlo simulation ClS db function of rj^ 

The first step uses the inclusive electron data to set the individual tower gains. Tower 
gains are determined in four time periods. The time boundaries correspond to natural breaks 
such as extended shutdowns or changes in accelerator conditions, so the statistics for each 
time period are not the same. The mean numbers of events per tower are 190, 190, 750, 
and 600, respectively, for the four time periods. These correspond to statistical precisions 
on the tower gain determination of ±0.64%, ±0.64%, ±0.33%, and ±0.38%, respectively. 

Having determined the individual tower gains, long-term drifts within each time period 
are measured by fitting to a line based on run number (typically a run lasts about 12 hours). 
These corrections remove aging effects or seasonal temperature variations, but are insensitive 
to short term variations such as thermal effects caused by an access to the detector in the 



The material traversed by electrons increases with polar angle, so {E/p) increases with \r] 
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collision hall. 

The next step uses the W sample to update the mapping corrections which describe 
the variation in response across the face of the towers. The strip chamber determines the 
local X (azimuthal) and z (polar) coordinates within the wedge, where —24 < x < 24 cm 
is measured from the tower center and —240 < z < 240 cm from the detector center. 
The (E/p) distribution as a function of x is fitted to a quadratic function, which corrects 
primarily for non-exponential attenuation in the scintillator of the light seen by the two 
phototubes. Tower-r^-dependent corrections are also made as a function of z. The statistical 
uncertainty in the mapping corrections is 0.2% in x and 0.13% in z. 

Finally a very small correction takes into account a systematic difference of the "under- 
lying event" in the inclusive electron and W datasets. The underlying event consists of two 
components - one due to additional interactions within the same beam crossing (multiple 
interactions) and the other due to the remnants of the protons and antiprotons that are in- 
volved in the inclusive or W electron production. It overlaps with the electron, contributing 
approximately 90 MeV on average to the electron Et- Because of the difference in Ej- be- 
tween the inclusive electrons (< Et >~ 10 GeV) and the W electrons (< Et >~ 38 GeV), 
their underlying energy contribution is proportionately different. This difference varies with 
the instantaneous luminosity, which is strongly correlated with time. 



All of the corrections applied to the W electrons are shown in Figure |T5|. The mean 
temporal correction is +4.6% and the mean mapping correction is —2.5%. The corrections 
reduce the RMS width of the E/p distribution from 0.0578 to 0.0497. 
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FIG. 15. (a) Spatial and temporal energy correction factors on the W electrons. The dotted 
curve shows the spatial corrections only, the dashed curve the temporal corrections only, and the 
solid curve the product of the two. (b) The E/p distributions of the W electrons after the respective 
corrections. The squares show the data before any corrections are applied. The improvement in 
the resolution after correction is apparent. 
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C. Event Selection: W eu, Z ^ e^e 

The W ^ eu and Z e~^e~ selection criteria are chosen to produce datasets with low 
background and well-measured electron energy and momentum. They are identical to those 
for the Z fi^fi^ and W ^ hp datasets except for the charged lepton identification and 
the criteria of removing Z e^e~ events from the W ^ eu candidate sample. The cuts 
and number of surviving events are shown in Table Vlll] and the electron criteria and the Z 



removal criteria are described in detail below. The samples begin with 108,455 W candidate 
events and 19,527 Z candidates events that pass one of two level-3 W oi Z triggers, and 
have an "uncorrected" electromagnetic cluster with Et > 20 GeV and an associated track 
with pt > 13 GeV/c. 

Candidate electrons are required to be in the fiducial region. This requirement primarily 
removes EM clusters which overlap with uninstrumented regions of the detector. To avoid 
azimuthal cracks, |x| is required to be less than 18 cm, and to avoid the crack between the 
2; > and z < halves of the detector, \z\ is required to be greater than 12 cm. The 
transverse EM energy is required to be greater than 25 GeV, and to have an associated 
track with pt > 15 GeV/c. The track must pass through all eight superlayers of the CTC, 
which improves the electron purity and limits the occurence of very hard bremsstrahlung. No 
other track with pj- > 1 GeV/c associated with the nominal vertex may point at the electron 
towers. This criterion reduces the QCD dijet background in the W sample. It also has the 
effect of removing the W and Z events which have secondary tracks associated with the decay 
electrons. These secondary tracks can result from the conversion of hard bremsstrahlung 
photons or through accidental overlap with tracks from the underlying event. Both of these 
sources are included in the simulation. Events are rejected when another track has an 
invariant mass below 1 GeV when combined with the electron cluster. 

A Z — > e^e^ event can fake a W ^ eu event if one of the electrons passes through a 
crack in the calorimeter. Most of these electrons are in the tracking volume. An event is 
considered to be a Z candidate if there is a second track with px > 10 GeV/c which has 



56 



opposite sign to the electron track and points at either the 6 = 90° or 6 = 30° crack, or is 
extrapolated to |a;| > 21 cm in the strip chamber. Z candidate events are removed from the 
W sample. For the Z sample, the two electron tracks are required to have opposite sign. The 



selection criteria described above are properly included in the Monte Carlo simulation pH 
The transverse mass in the region 65 < < 100 GeV/c^ and the invariant mass in 
the region 70 < M < 110 GeV/c^ are used for extracting the W mass and the Z mass, 
respectively. These transverse and invariant mass cuts apply only for mass fits and are 
absent when we otherwise refer to the W or Z sample. The final W sample contains 42,588 
events, of which 30,115 are in the region 65 < Mt < 100 GeV/c^. The final Z sample 
contains 1,652 events, of which 1,559 are in the region 70 < M < 110 GeV/c^. The E^, Elf, 



and Mt after all cuts are shown in Figure 16 for the W sample. 



Criterion 


W events after cut 


Z events after cut 


Initial sample 


108,455 


19,527 


Z vertex requirement 


101,103 


16,724 


Fiducial requirements 


74,475 


9,493 


Tracks through all CTC superlayers 


71,877 


8,613 


E^. > 25 GeV 


67,007 


6,687 


E!f > 25 GeV 


55,960 


N/A 


|u| < 20 GeV 


46,910 


N/A 


P| > 15 GeV 


45,962 


5,257 


-^tracks the electron towers = 1 


43,219 


1,670 


Me,track < 1 GeV 


43,198 


N/A 


Not a Z candidate 


42,588 


N/A 


Opposite sign 


N/A 


1,652 


Mass fit region 


30,115 


1,559 



TABLE VIII. Effect of selection cuts. 
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FIG. 16. Kinematic quantities from the final W ^ ev sample. Et distributions of (a) electrons 
and (b) neutrinos. The dashed curves show the events in 65 < Mt < 100 GeV, the fit region for 
the W mass measurement, (c) Transverse Mass distribution. The arrows indicate the region used 
in the W mass fit. 



D. Electron Energy Scale and Resolution 



All calibrations described above IV B are relative corrections designed to improve uni- 
formity. The energy scale is extracted from the reconstruction of the Z mass. The Z Monte 
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Carlo events are generated in the manner described in Section [III F| . The Monte Carlo events 
are then processed through the detector simulation where the electron energy is smeared 
according to the resolution: 



where all energies are in GeV, the stochastic term 13.5% was measured in the test beam, 
and the constant k includes such effects as shower leakage and residuals from the uniformity 
corrections discussed in Section 4.2. The parameter k is allowed to vary in the Z mass fit. 
The other variable parameter in fitting the Monte Carlo events to the data is a scale factor, 
Se- 

For the fit, a binned maximum likelihood technique is used where the data and Monte 
Carlo events for Mz are divided into 1 GeV/c^ bins for the interval 70 — 110 GeV/c^. The 
results are: 

Se{Z) = — 1^ = 1.0000 ± 0.0009 (10) 

Mz 

and 

K = (1.53 ±0.27)%, (11) 



where the uncertainties come from the Z statistics. The fit results are shown in Figure \17 . 
The two parameters are largely uncorrelated. The value of Se is equal to 1 by construction; 
the initial value of Se was not 1, but we iterated the fit with the scale factor applied to the 
energy until the final scale factor becomes 1. 



59 




M^^ (GeV/c^) 

^ 2.25 I 1 

S (b) 2o 

2 - 
1.75 
1.5 
1.25 
1 

I I I I I I I I I I I I I I I I I I I I I I I I 

0.997 0.998 0.999 1 1.001 1.002 1.003 

/ M/DG 

FIG. 17. Results of fit to Z mass and energy resolution, (a) Invariant mass distribution. The 
points are the data, and the solid line is the Monte Carlo simulation (normalized to the data) with 
best fit. (b) Correlation between the scale factor {Se) and the constant term (k) in the resolution 
function. 

A number of checks are performed to insure that these results are robust and unbiased. 
For example, 1000 Monte Carlo subsamples are created where each sample has the same 
size as the data, and are used to check that the likelihood procedure is unbiased and that 
statistical uncertainties by the fit are produced correctly. Moreover, compatible results 
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are found when a much simpler technique is used, that is, comparing the mean Mz, in 
the interval 86 — 96 GeV/c^, between the data and the Monte Carlo events. The Monte 
Carlo events include a 1% QCD background term. If the background term were omitted 
entirely, the energy scale and k would change by much less than their statistical uncertainties; 
we conclude that the uncertainties in the background have negligible contribution to the 
uncertainties in the fit results. Finally a Kolmogorov-Smirnov (KS) statistic is used to 
quantify how well the Monte Carlo events fit the data. The probability that a statistical 
fluctuation of the Monte Carlo parent distribution would produce a worse agreement than 
the data is 19%. The likelihood fit is also checked by varying the parameters in the KS fit 
to find a maximum probability. The result is Se = 1.0007 ±0.0010, in good agreement with 
the likelihood method. 



E. Energy Nonlinear ity Correction 

The average Et for Z decay electrons is about 4.5 GeV higher than those for W decay. 
Since the energy calibration is done with the Z's, any nonlinearity in the energy response 
would translate to an incorrect energy scale at the W. The nonlinearity over a small range 
of Et can be expressed as 

^ = e X AEr. (12) 

The slope, ^, could arise from several sources: energy loss in the material of the solendoid, 
scintillator response versus shower depth, or shower leakage into the hadronic part of the 
calorimeter. The near equality of the E/p scale factors for the W and Z samples limits the 
slope to be less than about 0.0004 GeV^^. The spread in electron Et for each of the W and 
Z samples is larger than the difference in the averages, so the most sensitive measure of ^ 
is the variation of the mean E/p between 0.9 and 1.1 for both samples as a function of Et- 



Their Et distributions and the residuals, (-E'/p)data— (-E/p) simulation, are shown in Figure [18 
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FIG. 18. Left: The Et distributions of electrons from W and Z decays. Right: Residual of 
data and Monte Carlo fit to E/p versus electron Et for the W and Z samples. The solid line is a 
linear fit with x^/dof = 1.4. When the slope is forced to be zero, the x^/dof increases to 2.2. The 
arrows represent the average Et values of the electrons for the W and Z samples. 

A linear fit to the E/p residuals for the W and Z data yields a slope of (1.91 ± 0.58) x 
10~^ GeV~^ in {E/p). Correcting the relationship between {E/p) and the scale factor 
gives a slope ^ = -0.00029 ± 0.00013(stat.) ± 0.00006(syst.) GeV-\ where the systematic 
uncertainty comes from backgrounds and the fitting procedure. The electron Et is corrected 
by 

Et Et{1 - 0.00029(Et - 42.73 GeV)) (13) 

before the final fit for the W mass. This correction shifts the fitted W mass up by (34 ± 
17) MeV/c^. The mean Et for the Z sample is 42.73 GeV, so the energy scale is unchanged 
at that point. 
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F. Check of Energy Scale and Momentum Resolution Using E/p 



The momentum scale was set with the Z mass as discussed in Section |T|, 

and the energy scale was set with the Z e+e~ mass as discussed in this section. In 
principle, the electron energy scale can be set by transferring the momentum scale from the 
T(ls) or J/ip ^ fi^fi~ mass as done in the Run lA analysis and equalizing E/p for data 
and simulation in W ^ ev decays. This technique has great statistical power and indeed 
was the preferred technique in previous CDF publications of the W mass @,|T^. However, 
systematic effects in tracking electrons are potentially much larger than for muons due to 
bremsstrahlung. To accurately simulate external bremsstrahlung effects ^Of' ^^e Monte 



Carlo program includes the magnitude and distribution of the material (see Section [TTTAD 



traversed by electrons from the interaction region through the tracking volume, propagation 
of the secondary electrons and photons,^^ and a procedure handling the bias on the beam 
constrained momentum which is introduced through the non-zero impact parameters of 
electrons that have undergone bremsstrahlung . 



To fit to the E/p distribution (see Figure |T^) to determine the energy scale, the width 



of the E/p distribution needs to be understood. It has a contribution from both the E 
resolution and the p resolution. At the W electron energies, the p resolution dominates. 
When the E/p distribution is fit to determine the energy scale, the E resolution is fixed to 
the value determined by the Z data, and the l/pr resolution is allowed to vary. As can be 
seen from Figure the E/p distribution agrees well with the resolution values determined 
solely from the Z yU+yU^ data. However, there is an excess at the low E/p tail region. 
Studies of the transverse mass for data events in this region show that the tail is due to 
mis-measured tracks in real W events. To account for this excess, the track parameters 
are smeared according to a second, wider Gaussian term for 8% of the Monte Carlo events. 
The two Gaussians describe the overall E/p distribution well. However, adding the second 



WXhe photons are treated in the same manner as the electrons in the calorimeter simulation. 



63 



Gaussian distribution does not significantly change the derived scale. 



10 W ^ev 




1 



0.8 0.85 0.9 



0.95 



1.05 



1.1 



1.15 1.2 

E/p 



100 







i 



■100 



0.8 0.85 0.9 



0.95 



1 



1.05 



1.1 



1.15 1.2 
E/p 

FIG. 19. Top: E/p distribution for W events (points) and the best Monte Carlo fit. The 
sohd histogram is the Monte Carlo fit normalized to data, and the points are the data. The fit 
reproduces the shape very well as indicated by the x^/dof =0.86. Bottom: The difference between 
the data and the best fit simulation. 
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FIG. 20. The energy resolution k and tracking resolution 6{1/pt) as determined from fits to 
the E/p distribution in ^ eu events, compared to the same resolutions determined from the 
Z ^ e^e~ and Z fJ-^f^~ data. 

The E/p distribution is fit for an energy scale and tracking resolution using a binned 
likehood method. The method is similar to the one used to fit the Z mass. The data are 
collected in 25 bins for the region 0.9 < E/p < 1.1, containing 22,112 events as shown in 



Figure |T9|. The log likelihood is maximized with respect to Se and the momentum resolution 
simultaneously. The energy scale factor is found to be 
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Se{E/p) = 0.99633 ± 0.00040(stat.) 

±0.00024(/t) ± 0.00035(Xo) ± 0.00018(pt scale), 

where 0.00024 comes from the uncertainty in the calorimeter resolution, 0.00035 from the 
uncertainty in the radiation length measurement, and 0.00018 comes from the uncertainty 
in the momentum scale which for this purpose is determined by the T(ls) measurement 
(see Section [III G|) . The result of the fit is shown in Figure 19. When we account for the 



nonlinearity of the calorimeter energy between Z decay electrons and W decay electrons as 
described in Section IV E| , the scale factor becomes 



Se{E/p) = 0.99480 ± 0.00040 (stat.) (14) 
± 0.00024 (k) ± 0.00035 (Xo) ± 0.00018 (pt scale) 
± 0.00075 (GEM nonlinearity). 

It is in poor agreement (3.9ct discrepant) with the energy scale determined from the Z mass 
(Eq. [Ty). When this scale factor is applied to the data, the Z mass is measured to be 0.52% 
lower than the world-average value. 

The E/p distribution for the Z sample is also used to extract Se- The result is: 

Se{E/p) = 0.99720 ± 0.00130 (stat.) (15) 
± 0.00024 (k) ± 0.00035 (Xo) ± 0.00018 (pt scale). 

The systematic uncertainties with respect to k, Xo, and momentum scale are common for 
the W and Z samples. The difference between this scale value and the scale from the Z 
mass is 2.0a. When both the W and Z events are combined, the discrepancy is 5.3a. 

The disagreement between the energy scale determined from the Z mass (Eq. |I0[) with 
that determined by the E/p distribution (Eq.s 14 and 15) is significant; therefore it would 
be incorrect to average the two. Moreover, the two techniques applied to the Z sample use 
the same energy measurements, thus hinting at a systematic problem between the tracking 
for muons and that for electrons, or a systematic difference between the actual tracking and 
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the tracking simulation. Another possibihty is an incomplete modeling of the calorimeter 
response to bremsstrahlung in the tracking volume. Appendix A describes some possible 
causes. 

As a result of this disagreement, we choose to use conservative methods for both the elec- 
tron energy and muon momentum scale determination. We use the Z — > e^e~ mass instead 
of the Ej]) distribution to set the electron energy scale since this is a direct calibration of 
the calorimeter measurement without reference to tracking or details of the bremsstrahlung 
process. Although statistically much less precise, we use the Z A*"*"/^" mass instead of the 
T(ls) or J 1'^ mass to set the muon momentum scale. 



The electron energy scale is determined by normalizing the measured Z ^ e^e" mass 
to the world-average mass. The measurement is limited by the finite statistics in the Z 
peak which gives the uncertainty of 72 MeV/c^ on M^r. A small nonlinearity is observed, 
resulting in /S.Mw — (34 ± 17) MeV/c^. Adding these uncertainties in quadrature, the 
total uncertainty on Mw due to the energy scale determination is 75 MeV/c^. The energy 
resolution is measured from the width of the Z — > e^e~ peak in the same dataset: — 



-I- (1.53 ± 0.27)%^. When the electron energy resolution is varied over this allowed 



Backgrounds in the W samples come from the following processes: 

1. W ^ TV ^ Euvu 

W —>■ Tu ^ hadrons -|- 

2. Z — > i'^i" where the second charged lepton is not detected 

3. Dijets (QCD) where jets mimic leptons 

4. cosmic rays 



G. Summary 




range in the simulation, the measured W mass changes by 25 MeV/c^. 



V. BACKGROUNDS 
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Contributions from Z ^ t^t , W^W , and tt are negligible. In general, backgrounds have 



a lower average transverse mass than W ^ iu decay, and, if not accounted for, will lower 



the fitted mass. All the background distributions as shown in Figure ^ are included in the 
simulation. 
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FIG. 21. The fractions (%) of backgrounds as a function of transverse mass distribution for the 
W ^ ev sample (Top) and the W ^ fiv sample (Bottom). The smallest contributor, W ^ tv ^ 
hadrons + vv, is not shown in this figure. 
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A. W ^ eu Backgrounds 

Few W TV ^ evvv events pass the kinematic cuts since the electron the total 
neutrino |Et|, and M-^ are substantially lower than those in the W ^ ev decay. W —* ru 
evvv events are estimated to be 0.8% of — > ev events in the W mass fitting region. This 
is the largest background in the W ^ ev sample, and is also the easiest to simulate. We 
have also simulated the W rv background where the r decays hadronically. We expect it 
to be (0.054 ± 0.005)% of the W sample. After Z removal cuts, very few Z —>■ e^e~ events 
can mimic W ^ ev events. The Monte Carlo simulation predicts (0.073 ± 0.011)% of the 
W sample in the mass fitting region to originate from Z e+e^. 

Dijet events can pass the W selection cuts if one of the jets mimics an electron and the 
other is mismeasured, creating ^x- Such events are refered to as "QCD" background. The 
QCD background is estimated by selecting QCD candidates from the W sample without 
Mj- and |u| cuts and plotting distributions of |u| and Mj- as shown in Figure (a detailed 



description can be found in Reference 1|28|). The number of QCD events predicted in the 
signal region "Region A" (see the top figure) is given by 

^j- _ Njiegion A (QCD) ,^ 

Region A (W) — T7 X Region B (W) 

Region B (QCD) 

= 249 ± 108, 

from which we find 119±56 events or (0.36±0.17)% of the W events are in the W mass fitting 
region. The kinematical distributions of the QCD events are derived from the W ^ ev 
sample with inverted electron quality cuts. 
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FIG. 22. Mt (GeV/c^) vs |u| distributions without Mr and |u| cuts for all W data (Top), and 
a QCD subset of the W data (Bottom). 



W ^ liv Backgrounds 

The largest background in the W ^ /lu sample comes from the Z /J'^/J'" process with 
one of the muons exiting at low polar angle (outside of the CTC volume) which mimics a 
neutrino in the calorimeters. The simulation predicts this background to be (3.6 ± 0.5)%. 
The uncertainty in the background estimate comes from two sources: the uncertainty in the 
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measured tracking efficiency at large 77, and the choice of parton distribution functions. 

The second largest background comes from the W tu process where r /iz/z/z/, 
which is 0.8% of the W sample. The W tv background where the r decays hadronically 
is negligible. Background from QCD is estimated by using the data in a similar manner 
to the electron case. The W ^ fiu sample is estimated to contain (0.4 ± 0.2) % of its 
events from the QCD process. Cosmic rays can appear as two oppositely charged back-to- 
back tracks in when they cross the detector in time with pp collisions. Most of them are 
removed by the W ^ selection criteria such as the Z removal cut or \Dq\ < 0.2 cm (see 
Section pi C|) . The number of cosmic rays remaining in the final sample is estimated by 
using events which fail \Dq\ < 0.2 cm criteria, but which pass all the other selection criteria. 
The expected number of cosmic ray events corresponds to (0.10 ± 0.05)% of the W sample. 



C. Summary 



Table |^ summarizes the fraction of the background events in the W samples in the 
mass fitting region. The total backgrounds in the W ^ ev and W —>■ fiu fit region are 
expected to be (1.29 ± 0.17)% and (4.90 ± 0.54)%, respectively. Adding the backgrounds in 
the simulation leads to shifts of (80 ± 5) MeV/c^ and (170 ± 25) MeV/c^ in the W ^ eu 
and W ^ fiu mass measurements, respectively. 



Background source 


W ^ ev sample 


W ^ fiv sample 


W ^ TV ^ ivvv 


0.8% 


0.8% 


W ^ TV ^ hadrons -|- vv 


(0.054 ± 0.005)% 




Lost Z 


(0.073 ±0.011)% 


(3.6 ± 0.5)% 


QCD 


(0.36 ±0.17)% 


(0.4 ±0.2)% 


Cosmic rays 




(0.10 ±0.05)% 


Total 


(1.29 ±0.17)% 


(4.90 ±0.54)% 



TABLE IX. Backgrounds in the W ^ ev and W ^ ^v sample in the mass fitting region. 
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VI. W PRODUCTION AND DECAY MODEL 



We use a Monte Carlo program to generate W events according to a relativistic Breit- 
Wigner distribution and a leading-order {p^ = 0) model of quark-antiquark annihilation. 
The distribution in momentum of the quarks is based on the MRS-R2 parton distribution 
functions (PDFs) The generated W is Lorentz-boosted, in the center-of-mass frame of 
the quark-antiquark pair, with a transverse momentum, pj^. The pj^ spectrum is derived 
from the Z e^e~ and Z fJ'^fJ'^ data and a theoretical prediction for the ratio of Z and 
W pt spectra which is differential in the rapidity of the vector boson. The Monte Carlo 
program also includes QED radiative effects . 

A. Parton Distribution Functions 

The uncertainty associated with PDFs is evaluated by varying the choice of PDF sets and 



by parametric modifications of PDFs. Figure 23 shows the CDF data on the W lepton charge 



asymmetry [^] which is sensitive to the ratio oidiou quark densities {d/u) at a given parton 
momentum fraction, x. Of all modern PDFs, the two giving the best agreement, MRST 



and CTEQ-5 |^ , are shown.^ Unfortunately the agreement even with these PDFs is barely 
satisfactory. Hence we follow reference [5D| in making parametric modifications to the MRS 



family of PDFs. These modifications with retuned parameters are listed in Table ^ and 
their predictions are compared to the W lepton charge asymmetry measurement and the 
NMC d/u data in Figure From the variation among the six reference PDFs, an 
uncertainty of 15 MeV/c^ is taken which is common to the electron and muon analyses. 



^^Predicted W charge asymmetries are calculated with the DYRAD NLO W production pro- 
gram isoi. 



72 




0.25 0.5 0.75 1 1.25 1.5 1.75 2 2.25 2.5 

I Lepton Rapidity I 



FIG. 23. The CDF W lepton charge asymmetry measurement compared to predictions using 
the DYRAD calculations with MRST (solid) and CTEQ-5 (dashed) PDFs. 



PDFs 


Modification 


MRST 

MRS-R2 

MRS-Rl 


d/u d/u X (1.07 - 0.07e-^^) 
d/u d/u + O.llx X (1 + x) 
d/n ^ d/n X (1.00 - 0.04e 0.015 )) 



TABLE X. Reference PDFs and modifications. 
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FIG. 24. (a) The CDF measurement of the W lepton charge asymmetry compared with the 
six reference PDFs. The upper and lower dotted curves are MRS-R2 and MRS-R2 modified, the 
upper and lower dashed curves are MRS-Rl modified and MRS-Rl, and the upper and lower solid 
curves in |?7| < 1 are MRS-T and MRS-T modified, respectively, (b) The NMC d/u data evolved 
to = M^. The gray bands represent the range spanned by the six reference PDFs. 
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B. W Transverse Momentum Spectrum 



The spectrum of W transverse momentum, , is needed to simulate the hneshape of 
transverse mass. The W mass measurement uses events at low pj^ where the theoretical 
calculations are not reliable. It would be difficult to extract pj^ from the W data because the 
neutrino momentum is not well measured. However one can model p^ through a measure- 
ment of , which can be measured accurately using the charged leptons from the Z decays. 
Theoretical calculations predict the cross-section ratio of W^s and Z's as a function of pt 
with small uncertainty since the production mechanisms are similar [Q. The measurement 
of is combined with the theoretical calculations of the ratio to derive p^ . This procedure 
is applied separately to the muon and electron samples, so the derived pj^ distributions are 
essentially independent although compatible. 

For each Z sample, a functional form for the Z px distribution is assumed for input to 
a Monte Carlo generator. The lepton response is modeled according to detector resolution 
and acceptance. The parameters of the assumed functions are fit to give agreement with 



the observed Z px distributions. The observed Z px distributions are shown in Figure ^ 
and are compared with the simulation which uses the best fit parameters for the input p^ 
distribution. 
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FIG. 25. The observed Z px distributions (points) for the (a) Z — > e^e~ and (b) Z — > fJ.~^n~ 
sample are compared with the Monte Carlo simulation. The solid line in (a) shows the best fit 
parameters for the input distribution, whereas the shaded band in (b) shows the la variation 
of the fit parameters. 



Resummed calculations are used for correcting the difference between the W and 

Z pt distributions, in terms of the ratio of the two distributions. As shown in Figure ^ (a), 
(b) and (c), the ratio is between 0.9 and 1.0 over the range of interest. Effects from the 
large ratio at ~ is very small since da/d{pT) a.s pr —* 0- The variation of the ratio 
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is studied by varying PDFs and nonperturbative parameters in the resummed calculations, 



and by calculating it in two different resummed schemes, one in impact parameter space |^ 
and the other in pt space . There is a rapidity dependence to the px distribution, 

illustrated in Figure ^ (d) and (e). This rapidity dependence is taken into account when 
pY is derived from p^. As indicated in Figure ^ the range of the possible ratio and rapidity 
dependence variation is about 2%. 

The extracted p^ distribution for the muon channel at the generation level is shown in 
Figure ^ (b). The shaded band represents the total uncertainty on the distribution. 
The dominant uncertainty comes from the finite statistics of the Z sample. The theoretical 
uncertainty in the pt ratio and rapidity dependence is small. The fractional uncertainties on 
the distribution from the statistics and theoretical calculations are shown in Figure ^ (a). 

The uncertainty on the W mass is evaluated by varying the p^ distribution within the 
shaded band in Figure (a). The finite statistics of the Z sample contributes independent 
uncertainties of 15 MeV/c^ and 20 MeV/c^ for the W ev and W ^ fiu channel. The 
contribution of the theoretical uncertainty is 3 MeV/c^ which is common for the electron 
and muon channel. 
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FIG. 26. The ratios of the pjf to distribution from resummed calculations in impact pa- 
rameter space showing (a) PDF dependence, and (b) nonperturbative parameter dependence. The 
ratios in impact parameter space and px space are compared in (c). The ratio of at = 0.3 
to for (d) 0<y^ < 0.25, and (e) 0.25 < < 0.5. 
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calculations, (b) The pj^ distribution extracted from the distribution and the theoretical 
calculations of Pj^ /Px for the W ^ jiv mass measurement. The band represents the uncertainties. 



C. QCD Higher Order Effects 



The W bosons are treated as spin-one particles and decay via the weak interaction into 
a charged lepton (e, ii or r) and a neutrino. The charged leptons are produced with an 

79 



angular distribution determined by the 0{al) calculation of which, for bosons with 
a helicity of -1 with respect to the proton direction, has the form : 

oc 1 + ai{pT) cos6'cs + 02 (pt) cos^ 9cs (16) 



d cos ^cs 

where px is the transverse momentum of the W and ^cs is the polar direction of the charged 
lepton with respect to the proton direction in the CoUins-Soper frame [|^. Oi and 02 are 



Pt dependent parameters. For px = 0, ai = 2 and 02 = 1 providing the angular distribution 
of a W boson fully polarized along the proton direction. For the pj^ values relevant to the 
W mass analysis {pj^ <~ 30), the change in W polarization as pj^ increases only causes 
a modest change in the angular distribution of the decay leptons [|^. The uncertainty is 
negligible. 

D. QED Radiative Effects 

W'j production and radiative W decays {W — > ih"y) are simulated using the calculation 



by Berends and Kleiss pT|j3^ . Most photons tend to be collinear with the lepton, often 
showering in the same calorimeter towers as the lepton. For the electron channel, these 
photons are merged with the electron cluster; for the muon channel, they reduce the muon 
momenta by their energy. Radiative effects from collinear photons are thus expected to 
be larger in the muon channel. Photons not collinear with the lepton are included in the 
calculation of u (see Figure 0), and have an effect that is similar in both the electron and 
muon channels. 

Shifts in the W mass due to radiative effects are estimated to be (—65 ±20) MeV/c^ and 
(—168 ± 10) MeV/c^ for the electron and muon channel, respectively. Uncertainties of the 
radiative effects are estimated from uncertainties in the theoretical calculation and in the 



calorimeter response to the photons. The Berends and Kleiss calculation [^] does not include 
all the radiative Feynman diagrams. For example, it does not include initial state radiation 
(if:- and w-channel diagrams) and allows a maximum of one photon. The effect arising 
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from the missing diagrams is evaluated by incorporating the PHOTOS package which 
allows two photon emissions, and the calculation by U. Baur et al. who have recently 
developed a complete 0{a) Monte Carlo which incorporates the initial state QED radiation 
from the quark-lines and the interference between the initial and final state radiation as well 
as including a correct treatment of the final state soft and virtual photonic corrections. The 
effects on My/ from the former case are less than 10 MeV/c^ for the W ^ ev channel and 
less than 5 MeV/c^ for the W ^ fiu channel. The effects on Mw from the latter case are 
less than 20 MeV/c^ for the W ^ ev channel and ~10 MeV/c^ for the W ^ fiu channel. 
The uncertainty in the calorimeter response to the photons well-separated from the W decay 
lepton, is evaluated by varying the photon energy threshold, the photon fiducial region, and 
the photon energy resolution. The effect is 3 MeV/c^ on the W mass. 



E. Summary 

The uncertainty associated with PDFs is evaluated by varying the choice of PDF sets. It 
is estimated to be 15 MeV/c^ which is common to the electron and muon analyses. The 
spectrum is derived from the Z — * e^e~ and Z fi^fi^ data and a theoretical prediction 
for the ratio of Z and W pr spectra differential in the rapidity of the vector boson. The 
corresponding uncertainty in the W mass is dominated by Z statistics. It is 15 MeV/c^ for 
the W ^ ev channel and 20 MeV/c^ for the W —>■ fiu channel. A common uncertainty of 
3 MeV/c^ comes from the theoretical prediction for the ratio. The uncertainty in the W 
mass due to QED radiative effects is estimated to be 20 MeV/c^ to the W ^ eu channel, 
and 10 MeV/c^ to the W fiu channel. 



VII. RECOIL MEASUREMENT AND MODEL 

The transverse mass distribution used for the W mass measurement is reconstructed 
using the Ex of the charged leptons (described in Section Q and fV\) and the neutrinos. 
The transverse energy of the neutrino is inferred from the charged lepton Et and the recoil 
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energy u (see Figure 0). This section describes the reconstruction of u, and an empirical 
model of the detector response to u which is implemented in the simulation. Since the W 
and Z share a common production mechanism and are close in mass, the recoil model is 
based mainly on Z decays. 



A. Recoil Reconstruction 

The recoil vector u is calculated by summing over electromagnetic and hadronic calorime- 
ter towers within the detector range |?7| < 3.6, 



u = (m^, Uy) = Etowers-E" siu 6'(cos (f), siu , 



(17) 



Table ^ lists tower thresholds for online (Level-3) reconstruction and this analysis. The 
thresholds for this analysis correspond to 5 times the calorimeter noise level. 



Calorimeter 


Online threshold (GeV) 


Analysis threshold (GeV) 


Central EM 


0.1 


0.1 


Central Had. 


0.1 


0.185 


Plug EM 


0.3 


0.15 


Plug Had. 


0.5 


0.445 


Forward EM 


0.5 


0.2 


Forward Had. 


0.8 


0.73 



TABLE XL Tower energy thresholds used to reconstruct u both in online and in this analysis. 
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There are two contributions to the recoil vector u. The first contribution is the energy 
of the initial state gluons radiated from the quarks that produce the W or Z boson. This 
energy balances the pt of the boson. The second is the energy associated with multiple 
interactions and the remnants of the protons and antiprotons that are involved in the W or 
Z production. The latter energy is referred to as the underlying energy. It is manifested in 
T^Et, where 

HEt = Stowers-S sin ^ = Etowers-^T- (18) 

The lepton energy should not be included in the u calculation, and thus the towers 
containing energy deposited by the lepton are excluded in the sum. This procedure removes 
two towers for muons, and two or three towers for electrons. If the center of the electron 
shower is more than 10 cm away from the azimuthal center of the tower {\x\ > 10 cm), there 
will be leakage in the azimuthally adjacent towers which are also removed. This procedure 
removes not only the lepton energy, but also the underlying energy which needs to be added 
back to the sum. The underlying energy is estimated from the energy in calorimeter towers 
away from the lepton in the W data. In the muon analysis, this energy is added back to the 
u calculation. In the electron analysis, rather than correcting u, the same amount of energy 
is removed from the Monte Carlo simulation. 

B. Recoil Model 

For the purposes of modeling the response and resolution, it is natural to define u in 
terms of the components Ui and U2, anti-parallel and perpendicular to the boson direction, 
respectively. The average value of ui is the average calorimeter response balancing the boson 
Pt, and the average value of U2 is expected to be zero. Ui and U2 are parameterized in the 
form 



/ \ 

Ui 






+ 


^^1(^71) 


\ 


[U2) 




) 




^ G2{a2) 


J 



(19) 
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where Gi{ai) and 6*2 (o'2) are Gaussian distributed random variables of mean zero and widths 
0"! and o"2, and the quadratic function is the response function to the recoil energy. 

A detailed description can be found in Reference |^ . 

The resolutions o"i and a2 are expected to be dependent on HEt- For the minimum bias 
events which represent the underlying event in the W and Z sample, the resolutions (ax) 



and (cXy) are well parameterized with TiEt- A fit to the data, as shown in Figure gives 

ambsi^Er) = 0.324 x (SEr)"-'"" (20) 

where ambsi'^Ex) and T,Et are calculated in GeV. For the W and Z events, a good de- 
scription of the resolution requires additional parameters which accont for its boson 
dependence; the initial state gluons balancing the boson px produce jets which contribute 
to the resolution differently than the underlying energy. In order to allow this resolution 
difference, the widths are parameterized in the form 



for the electron channel and 

(.A 

0"2 



(^mbsi'^Ex) X 



1 + Si ■ 

1 + S2- {p^r^'°''f 



(21) 



ambsi^Ex) X 



yo2y 



boson \ 



«i + A ■ p't' 

a2 + (32- P^t""" 



(22) 



for the muon channel, where Si, S2, ai, 02, Pi, and P2 are free parameters. Although the 
two channels use different formulae, the fitted funtions are consistent with each other - ai 
and 0^2 are close to 1 and the difference between the linear term and the quadratic term is 
within the statistical uncertainty of the Z sample. The argument J^Et in Eq.s ^ and ^ 
comes from the T^Erp distributions of the W and Z data. The T^Erp distributions in various 
bins are shown in Figure EUl They are nicely fit to F-distributions 



7(SET;a,6) 



\b-l^-a{EET) 



m 



(23) 



where a and b are fit parameters, and 6 is a linear function of The term a/T{b) 

normalizes the distribution. Figure ISO shows the J^Et distributions and fits for the Z and 
W events. 
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FIG. 28. The fit for the rms of the Ux and Uy distributions as a function of HEt using the 



minimum bias sample. 
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FIG. 29. The HEt distributions in 5 different bins for the Z ^^fi data are shown: (a) 
for <5 GeV, (b) for 5 < p| < 10 GeV, (c) for 10 < p| < 20 GeV, (d) for 20 < p| < 30 GeV, 
and (e) for 30 < p| < 50 GeV. 
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FIG. 30. The T,Et distributions for (a) the W ^ ev sample, (b) the Z — > e^e sample, 
the W ^ fiv sample, and (d) the Z IJ-^ fJ-^ sample. The solid lines are fits to the funct' 



described in Eq. 23 
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FIG. 31. The (ui) versus (solid lines) as derived from Z sample fits for (a) the electron 
channel and (b) the muon channel. The fits are compared with the data points. 
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FIG. 32. The fit of tlie ui and U2 distributions in ranges of Z pT in the muon sample, illustrating 
the adequacy of assuming Gaussian resolution(solid lines). 



The Z data provide ui, M2, SE't, and the pr oi the Z. The parameters in Eqs. 23, 21 



and |2^ are derived by fitting to these variables. Figure ^ compares {ui) as a function of 
from the Z data with the fit functions /{pt) described in Eq. |19[ The validity of a 



Gaussian parameterization in Eq. |T9| is illustrated in Figure The parameterization of the 
recoil response model is further cross-checked by distributions of ui, U2, and |u|. As shown 
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in Figure they all agree well. The u resolutions in the Z fi~^fi~ data are shown as a 
function of in Figure 0, where the data is compared with the recoil model with (the solid 
histograms) and without (the dashed histograms) including the effect of gluons against the 
W. As expected, the resolution gets worse in ui as the jet structure of the recoil becomes 
apparent, increasing T^Ej- in the Ui direction. 

While the Z sample, where the boson pt is well understood, allows the unfolding of 
response and resolution, the W samples do not allow these effects to be separately under- 
stood. However, the W samples can be used to optimize the model parameters for the W 



data while preserving a good description of the Z data. This is demonstrated in Figure |35 . 
The ultimate recoil model includes the |u| and u^ (the component of u perpendicular to the 
lepton direction) distributions from the W data in the fit. 
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FIG. 33. (a) ui, (b) U2, and (c) |u| distributions for the Z — /x+/x data. The histograms are 
the simulation using the recoil model parameters. 
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FIG. 34. a{ui) and (7(^2) as a function of pf for the Z sample. The points are the 

data, and the sohd histograms are the simulation using the recoil model parameters. The dashed 
histograms show ambsi^Ex), the resolutions of the underlying energy. 
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FIG. 35. The muon Z fits separately constrain resolution and response, as shown by the ellipse, 
while the W data gives a further correlated constraint, as shown by the band. This is obtained 
from the Monte Carlo studies. 

C. Comparison of Data and Simulation in the W Samples 

This section compares the data with the simulation which uses the best fit parameters 
of the modeUng. The W data is more naturally described in terms of components uy and 
u^ of recoil defined with respect to the charged lepton direction - the component along 
the lepton direction and the component perpendicular to the lepton direction, respectively 
(see Figure [36|) .p**! The |u| and uii distributions and residuals are shown in Figure |33 and 



Figure |3g. The u^ distribution is shown in Figure ^ The means for u^ are consistent with 



zero and the other u projection numbers are listed in Table The models reproduce the 
basic characteristics well. 



**When |u| << Ej^, the transverse mass becomes ~ 2Ej^ + uii. 



93 




V 



FIG. 36. Kinematics of leptons from the W decay and the transverse energy vector recoihng 
against the PF, as viewed in the plane transverse to the antiproton-proton beams, uy is the 
component of u along the lepton direction and the component of u perpendicular to the lepton 
direction. 
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FIG. 37. The (a) |u| and (c) uy distribution distribution for the W —i- eu sample. The points 
(histograms) are the data (simulation). The differences between the data and the simulation are 
shown in (b) and (d). 
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FIG. 38. The (a) |u| and (c) uy distribution for the W ^ jiu sample. The points (histograms) 
are the data (simulation). The differences between the data and the simulation normalized by the 
statistical uncertainty are shown in (b) and (d). 
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FIG. 39. (a) The distribution for the W ^ /iv sample. The points (histogram) are the data 
(simulation) . (b) The difference between the data and the simulation normalized by the statistical 
uncertainty. 



Quantity 


Mode 


Data 


Simulation 




eu 


5.684 ± 0.034 GeV 


5.765 GeV 






5.640 ± 0.065 GeV 


5.672 GeV 


^™(u||) 


eu 


5.877 ± 0.024 GeV 


5.827 GeV 


(^""""(uil) 




5.732 ± 0.069 GeV 


5.750 GeV 


(-11) 


eu 


-0.573 ± 0.034 GeV 


-0.639 GeV 


(-11) 




-0.436 ± 0.048 GeV 


-0.422 GeV 



TABLE XII. Widths and means for recoil response projections for data and simulation. The 
simulation includes the W constraint and background bias. Uncertainties shown here are only 
statistical, and do not include systematic uncertainties due to pj!^ and the recoil model. 
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One can further examine whether or not the model describes correlations among vari- 
ables. The distributions in un are examined in four bins of |u|, shown for the electron analysis 



in Figure ^ and for the muon analysis in Figure The correlation of uy and transverse 
mass is illustrated in Figure ^ and the trend of (uy) with azimuthal angle between the lep- 
ton and u is shown in Figure ^ As indicated in these figures, the simulation well represents 
the data. 
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FIG. 41. The uii distributions for the W ^ jiu sample in four bins of |u|. The points are the 



data, and the histograms the simulation. 
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FIG. 42. (a) The average value of uy as a function of Mt for the W ^ eu sample. The points 
are the data, and the solid histogram is for the simulation, (b) Residuals between the data and 
the simulation. 
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FIG. 43. (u||) as a function of azimuthal angle between the lepton and u for the W ^ fiv 
sample. The points are the data and the histogram is the simulation. 

D. Uncertanties on Mw 

The uncertainty on the W mass is evaluated by varying the model parameters within 
their uncertainties. The size of the parameter uncertainties is taken from the Z statistics 
and does not include the reduction produced by including the W data in the model. For 
each set of model parameters a set of transverse mass templates are produced which are fit 
to the transverse mass distributions of the data and a standard Monte Carlo template. The 
rms of Mw values obtained from the fit to the Monte Carlo template is 37 MeV/c^ for the 
electron channel and 35 MeV/c^ for the muon channel. 

E. Summary 

The detector response to the recoil energy against the W is modeled primarily using 
the Z data. The W data are used to optimize the model. The model is empirical 

in the sense that its form is justified by the data and its parameters determined from the 
data. The modeling procedure is applied separately to the muon and electron samples, so 
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the uncertainties on the W mass due to the recoil model are essentially independent. The 
parametrizations are compatible in the two channels. 

The uncertainty on the W mass is evaluated by producing a set of transverse mass tem- 
plates with the model parameters allowed within their uncertainties, and fitting to the trans- 
verse mass distributions of the data and a standard Monte Carlo template. It is 37 MeV/c^ 
for the electron channel and 35 MeV/c^ for the muon channel. 



VIII. RESULTS AND CONCLUSIONS 

This section summarizes the W mass results. Cross-checks which support the results 
are discussed. The results of the two lepton channels are combined with previous CDF 
measurements. The combined result is compared with other measurements and with global 
fits to all precise electroweak measurements which predict a W mass as a function of the 
Higgs boson mass. 



A. Fitting Procedure 

The W mass is obtained from a binned maximum likelihood fit to the transverse mass 
spectrum. This spectrum cannot be predicted analytically and must be simulated using 
a Monte Carlo program which produces the shape of the transverse mass distribution as 
a function of M^. This program incorporates all the experimental effects relevant to the 



analysis, including W production and decay mechanisms as described in Section VI, the 



detector acceptance for the charged leptons from the W decay, the detector responses and 



resolutions of the leptons as described in Sections 111 and K^, and the detector response 



and resolution of the recoil energy against the W as described in Section |V11[. The Monte 



Carlo program generates Mt distributions used as templates for discrete values of M^. The 
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width of the W, Tw, is taken as the Standard Model value |^ for that W mass.p^ The 
transverse mass distribution templates also include the background contributions. The mass 
fit compares the data transverse mass distribution to the templates. 

The transverse mass fitting procedure is tested by using large Monte Carlo samples and 
by generating pseudo-samples of the size of the data and extracting a mass value for each 
dataset. We investigated the bias in the fit and confirmed the statistical errors returned by 



the fits. The results are illustrated for the muon fit in Figure No biases are observed in 
the fitting procedure and the fit errors returned by the simulation datasets and the variation 
in returned mass values are consistent with the statistical uncertainties of the fits to the 
data. 



^I^^r^^ is precisely predicted in terms of the masses and coupling strengths of the gauge bosons. 
The leptonic partial width T{W iv) can be expressed as Gi? M^/6\/27r(l + 5sm) where 6sm is 
the radiative correction to the Born-level calculation. Dividing the partial width by the branching 
ratio, Br{W iu) = 1/(3 + 6(1 + as{Mw)/T^ + ©(a^))), gives the SM prediction for Tw 
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FIG. 44. (a) Difference between the input Mw values and the returned values by fits to Monte 
Carlo pseudo-samples. Each sample is 100 times the size of the W ^ data, (b) The (statistical) 
error returned by fitting 1000 Monte Carlo pseudo datasets of the same size as the W —>■ fiv data. 
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B. The W Mass Measurement 



The fit results yield the measurements of the W mass in the electron and muon channels. 
They are: 

= 80.473 ± 0.065 (stat.) ± 0.092 (syst.) GeV/c^ 



and 



Mi;^, = 80.465 ± 0.100 (stat.) ± 0.103 (syst.) GeV/cl 



as a 



The negative log likelihood distribution for the muon sample is shown in Figure ^ 
function of M^. A similar distribution is obtained for the electron sample. The transverse 
mass distributions for the W ^ ev and W ^ fiu samples are compared to the simulation 
with the best fits in Figures ^ and The fit curves give x^/dof of 32.4/35 and 60.6/70 
for the electron and muon samples, respectively. If we extend the region of comparison from 
65 < Mt < 100 GeV/c2 to 50 < Mt < 120 GeV/c^, the curves give xV^of of 82.6/70 and 
147/131, and Kolmogornov-Smirnov (KS) probabilities of 16% and 21%. 
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FIG. 45. The deviation of the negative log hkelihood from the minimum for the W ^ jiu 
sample. The W width is fixed at the Standard Model value in the fit. 
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FIG. 46. W Transverse mass distributions compared to the best fit for tfie W ^ eu cliannel. 
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FIG. 47. W Transverse mass distributions compared to the best fit for tfie W ^ fiv cliannel. 

A summary of all systematic uncertainties is given in Table 8.1. They are estimated by 
measuring the subsequent shifts in Mw when each source is varied by its uncertainty in the 
Monte Carlo simulation. The largest uncertainties come from the finite statistics of the Z 
samples. The Z statistics are the predominant source of the uncertainties on lepton scale, 
lepton resolution, the model, as well as the recoil model. As muon and electron analy- 
ses use the muon and electron Z sample separately, the statistical effects are independent. 



108 



The theoretical uncertainty in the pj^ distribution gives a small common-contribution. The 
uncertainty due to the choice of PDFs is evaluated for the muon acceptance and is essen- 
tailly the same for the electron acceptance. We take the PDF uncertainties to be identical 
and common for the two channels. Although the QED corrections are rather different for 
electrons and muons, there is common as well as independent uncertainty. 

The total common uncertainty for the two lepton channels is 16 MeV/c^, due almost 
entirely to the common determination of the parton distribution function contribution. Ac- 
counting for the correlations, the combined value is: 

Mw = 80.470 ± 0.089 GeV/c^ 



Source of uncertainty 


W ^ ev 


W ^ Hv 


common 


Lepton scale 


75 


85 




Lepton resolution 


25 


20 




PDFs 


15 


15 


15 


jdW 

r-j, 


15 


20 


3 


Recoil 


37 


35 




Higher order QED 


20 


10 


5 


Trigger & Lepton ID bias 




15010 




Backgrounds 


5 


25 




Total 


92 


103 


16 



TABLE XIIL Systematic uncertainties in the W mass measurement in MeV/c^. 
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FIG. 48. The 1-a and 2-a contours in Tw versus M]y of the transverse mass fit when the width 
is floated for (a) the W ev channel and (b) the W channel. The dashed lines are the 

predicted Tw as a function of Mw- 
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C. Cross-Checks of the W Mass Measurement 



The reliability of the measurement can be checked by fitting lepton pt instead of trans- 
verse mass, by subdividing the W samples, and by removing the constraint on the W width 
function of mass. 

The W width, T^, can be extracted from the transverse mass distributions by fitting 
either in the region near the Jacobean edge or in the high- Mr region. The CDF experiment 
measured Tw to be 2.04 ±0.14 GeV using 100 < Mt < 200 GeV/c^ By generating Mt 
templates at discrete values of and T^, and allowing them to vary in the fit, one can 
measure both Mw and Tw simultaneously from the region near the Jacobean edge. Since 
Tw provides similar effects to the input ^^^1 the detector resolution of u in this region, 
the measurement of Tw provides a check on the recoil and pj^ models. Figure ^ shows 
the 1-0" and 2-a contours of the fitted W width versus W mass. The widths are consistent 



with the Standard Model: it is almost identical to the SM value for the muon channel, and 
about 1.5 cr away for the electron channel. The fitted W mass differs by 60 MeV/c^ for 
the electron channel and 10 MeV/c^ for the muon channel from the values with F^^ fixed. 
We do not derive measurements of the width from these fits due to the large systematics 
variations which come from changing resolutions and modeling. 

The transverse momentum spectra of the leptons as shown in Figures and BO also 
contain W mass information. W mass values obtained from maximum likelihood fits are 
consistent with the values from the transverse mass fit. The distributions from the simulation 
with the best fits are compared with the data in the figures. 

The W mass results are cross-checked by making various selection criteria on the data and 
Monte Carlo simulation, and refitting for the W mass. The events are divided into positively 
and negatively charged lepton samples. For the electon sample the charge difference listed 
in Table pCIV| involves statistical uncertainty only and corrreponds to the mass difference of 
123 ± 130 MeV/c^ between the and the W~. For the muon sample the table entries 
include the tracking alignment uncertainty of 50 MeV/c^. The mass difference of 136 ± 205 



111 



MeV/c^ is observed between the and the W . The electron and muon results are 
combined to give a mass difference of 127± 110 MeV/c^. 



The samples are also partitioned into four bins of |u| as shown in Figures ^ and The 
Monte Carlo simulation reproduces the data very well in all the |u| bins, indicating that the 
W pt and recoil energy are well modeled in the simulation. When the events are partitioned 
into > 35 GeV/c and p^ < 35 GeV/c samples, the Mt shapes between the two samples 
(see Figure ^3]) are dramatically different. Yet there is good agreement between the data 
and simulation. 



The extracted W masses described above are summarized in Table XIV. 
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FIG. 49. Et distributions of (a) electrons and (c) neutrinos in the W ^ ev channel. The 
points are the data and the histograms the best fit simulation. The differences between the data 
and simulation are shown in (b) and (d). 
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FIG. 50. (a) pt distribution of muons and (c) Et distribution of neutrinos in the W iiv 
channel. The points are the data and the histograms the best fit simulation, (b) and (d) The 
difference between the data and simulation normalized by the statistical uncertainty. 
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Uj (GeV) (1 0<U<1 5 GeV) (GeV) (1 5<U<20 GeV) 

FIG. 51. Transverse mass distributions in bins of |u| for the W ^ eu data (triangles) and 
the best fit simulation (histograms). The four |u| bins are, < |u| < 5 GeV (Top Left), 
5 < |u| < 10 GeV (Top Right), 10 < |u| < 15 GeV (Bottom Left), and 15 < |u| < 20 GeV 
(Bottom Right) 
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FIG. 52. Transverse mass distributions in bins of |u| for the W ^ fiu data (points) and the best 
fit simulation (histograms). The four |u| bins are, < |u| < 5 GeV (Top Left), 5 < |u| < 10 GeV 
(Top Right), 10 < |u| < 15 GeV (Bottom Left), and 15 < |u| < 20 GeV (Bottom Right) 
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FIG. 53. Transverse mass distributions for (a) low pr and (b) high pT muons in the W ^ ^iv 
data (squares) and simulation (lines). 



Fitting 


Selection 


AMw (MeV/c^) (eu) 


AMw (MeV/c2) {/^u) 


Pj, 




-80 ± 60 


-19+132 


H/rp 




+76 ± 60 


-20 + 127 


Mt 




+62 ± 90 


+67 ± 145 
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Mt 


£- 


-61 ± 90 


-69 ± 145 


Mt 


< |u| < 5 GeV 


-1±86 


-41 ± 135 


Mt 


5 < |u| < 10 GeV 


-36 ± 110 


-164 ± 169 


Mt 


10 < |u| < 15 GeV 


+161 ± 204 


+484 ± 301 


Mt 


15 < |u| < 20 GeV 


-348 ± 385 


+534 ± 450 



TABLE XIV. Difference from the nominal value of extracted M]y values from lepton transverse 
momentum fits and from various subsample transverse mass fits. 



D. Combined W Mass 

The issue of combining the present results with previous CDF measurements 0] merits 
some additional discussion since the lepton energy and momentum scales were determined 
differently. In particular, in our the previous analyses the electron scale was determined 
with the E/p method. In the present work that procedure is shown to result in a Z mass 
discrepant by (0.52 ± 0.13)%; in the Run lA analysis, the discrepancy was (0.28 ± 0.24)%. 
The statistics of Run lA are insufficient to distinguish the two cases - that the E/p method 
worked well or was systematically off as indicated in the Run IB result. Moreover, the 
experimental conditions differ for the two runs. For example, the aging and rate effects in 
the CTC due to higher luminosity are more pronounced for the present work. For these 
reasons and because the underlying cause for the E/p discrepancy remains unresolved, we 
believe that applying a correction factor to the Run lA result is not warranted. We prefer 
to average the results as published with the stated errors. Thus the combined CDF result 
is: 

Mw = 80.433 ± 0.079 GeV/c^ 
This value is precise to 0.1% and corresponds to a total integrated luminosity of ~105 pb^^. 
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E. Comparison with Other Results 



The present results are compared with other pubhshed results in Table pCV| [^,^[^,0. 
The agreement is excellent. The direct measurement of the W mass is an important test 
of the Standard Model. The W mass is indirectly predicted precisely by including loop 
corrections involving the top quark and Higgs boson. The corresponding implication for the 
Higgs boson mass is shown in Figure Our result agrees well with the Standard Model, 
and when combined with all other electroweak results H] prefers a light Higgs boson. 



UA2 


80.360 ± 0.370 GeV/c^ 


CDF 


80.433 ± 0.079 GeV/c^ 






DO 


80.474 ± 0.093 GeV/c^ 






ALEPH 


80.418 ± 0.076 GeV/c^ upto . 




189 GeV 




80.423 ± 0.123 GeV/c^ upto , 


A = 


183 GeV 


DELPHI 


80.270 ± 0.144 GeV/c2 upto , 




183 GeV 


L3 


80.610 ± 0.150 GeV/c2 upto , 


A = 


183 GeV 


OPAL 


80.432 ± 0.080 GeV/c^ upto . 


/~s = 


189 GeV (preliminary) 




80.380 ± 0.130 GeV/c^ upto . 


/~s = 


183 GeV 


Indirect Meas. 


80.381 ± 0.026 GeV/c^ 



TABLE XV. Measurements of the W mass. CDF and DO measurements have a common error 
mostly due to Parton Distribution Functions. The LEP II measurements have common errors 
including the LEP beam energy. The indirect measurement includes the LEP and SLC Z pole 
measurements, the z^N measurement, and the Tevatron Top mass measurements. 
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F. Conclusions 



We have measured the W mass to be M^y = 80.470 ± 0.089 GeV/c^ using data with 
an integrated luminosity of ~85 pb^^ collected from 1994 to 1995. When combined with 
previously published CDF data, we obtain Mw = 80.433 ± 0.079 GeV/c^. 
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FIG. 54. The direct measurements of the W and top quark mass from CDF and D0 experi- 
ments, the direct measurement of the W mass from LEP II experiments, and the indirect W and 
top mass measurement from LEP, SLC, and Tevatron neutrino experiments. The curves are from a 
calculation of the dependence of the W mass on the top mass in the Standard Model using several 
Higgs boson masses. The band on each curve is the uncertainty obtained by folding in quadra- 
ture uncertainties on a^M^), Mz, and as{M'^). The uncertainty is dominated by the hadronic 
contribution to a(M|), Aahad = 0-028 ± 0.0007 (Ref.48). 
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APPENDIX A: DISCUSSION OF DISCREPANCY BETWEEN Mz AND E/P 

METHODS 

The calorimeter energy scale for the W mass measurement in this paper is set using the 
invariant mass distribution of Z — > e^e~ events. Ideally, the E/p distribution would be 
used to set the energy scale where the momentum scale is determined by the T — > ix^ix" 
data. The E/p distribution has a smaller statistical uncertainty than the method of using 
the Z — > e+e~ mass because it makes use of the higher statistics of the W and T samples. 
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The E /p method, however, gives a significantly different result than the Z — >■ e^e~ mass 
method. 

The Z — > e^e~ mass method gives the energy scale of 1 by construction (see Sec- 
tion fvg) : 

Se = = 1.0000 ± 0.0009. 

The E/p distribution for the W ^ ev data does not agree with the simulation with the 
energy scale given by the Z mass method. The best fit between the data and the simulation 
requires an energy scale, 

Se = 0.99613 ± 0.00040 (stat.). 



Including the non-linearity correction described in Section IV E the energy scale becomes 



Se = 0.9948 ± 0.00040 (stat.) 

± 0.00024 (k) 

± 0.00035 (Xo) 

± 0.00018 (pt scale) 

± 0.00075 (CEM Non - linearity) 

where the uncertainty on the momentum scale comes from the T mass measurement (see 



Section [III G\ ). The difference between the Mz result and the E/p result is 

1.0000 - 0.9948 



3.9 (Al) 



VO.OOO92 + 0.00102 

standard deviations. This is unlikely to be a statistical fluctuation. A Kolmogorov-Smirnov 
statistic is calculated for the comparison of the data to the Monte Carlo. The probability 
that a statistical fluctuation would produce a worse agreement in the integrated distributions 
is 5.5 X 10-^. 

This Appendix discusses checks given by various data samples, and possible explanations 
of the discrepancy between E/p and Mz methods. 
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1. Checks on E and p Scales 



The energy scale, Se, is checked using various data samples. The Z e~^e~ sample 
is used for extracting the E scale from E/p. The J/ip ^ and Z — ^ samples 

are used for extracting the p scale. The momenta of electron tracks for the ip e^e~ , 
T e+e~, and Z e^e~ samples are used for setting the p scale (see Figure The 
results are summarized in Table X.V1 and Figure While all the results are consistent with 
each other, the central values are closer to 1 when the E/p scale is determined using the 
Z ^ e^e~ sample instead of the W ^ ev sample, or when the p scale is determined using 
electron tracks instead of muon tracks. Problems in the electron non-linearity correction or 
differences between the electron and muon tracks beyond our simulation could cause this. 
However our results are not statistically significant enough to be conclusive. 
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FIG. 55. Invariant mass distributions of electrons using their momenta for V e+e~, 
T — > e+e", and Z — > e+e" data samples. The solid lines are the best fits from the Monte Carlo 
simulation. 



# 


Data Sample for p scale 


Data Sample for E/p scale 


Se 


Dev. from 1 
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1 


T — > LL^ LL 


W ^ eu 


0.9948 ± 0.0010 ± 0.0002 


-3.9a* 


2 


T — > UL^ u~ 


e+e~ 


0.9972 ± 0.0014 ± 0.0002 


-2.0(7 


3 


J lib — > 11^ ii~ 


W ev 


0.9947 ± 0.0010 ± 0.0004 


-3.8<T* 


4 


Z fi+fi- 


W ^ eu 


0.9952 ± 0.0010 ± 0.0011 


-2.8(T* 


5 


Z — > e"^e~ (tracks) 


W ^ ev 


0.9955 ± 0.0010 ± 0.0026 


-1.5(7* 


6 


T e"'"e~ (tracks) 


W ^ ev 


0.9970 ± 0.0010 ± 0.0020 


-1.2(7* 


7 


J/il^ — e+e~ (tracks) 


W ^ eu 


0.9959 ± 0.0010 ± 0.0015 


-2.0(7* 



TABLE XVI. Required energy scales for various data samples. The errors on Se come from 
the E/p scale (first) and the p scale (second). *: the deviation from 1 includes the Z statistical 
uncertainty (±0.0009). 
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FIG. 56. Required energy scales for various data samples. The shaded area represents the 
energy scale determined by the Z — > e~^e~ mass. 

2. Momentum Non-Linearity 



A non-linearity in the pt measurement could produce a discrepancy between the two 
methods. The average px of T (ip) decay muons is ~ 5.0 GeV/c (~ 3.5 GeV/c), while 
the average pt of W and Z decay electrons is ~ 40 GeV/c. Figure shows the difference 
between the measured mass and the expected mass as a function of the sum of l/pr of 
the two muons in T and ip decays. W and Z events occur on the far left of the plot. No 
significant momentum non-linearity is observed.^ 



ttt Without the new CTC calibration and alignment for this analysis, there appears to be a small 
non-linearity in momentum measurement (0.1% non-linearity from 2 GeV to 50 GeV). This went 
away with the CTC calibration and alignment. The change has not been fully understood. 
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3. Differences between the Electron and Muon Tracks 



In the E/p method, the electron momentum scale is determined from the muon momenta. 
In many ways, electron tracks are different from those of muons. They are produced with 
different internal bremsstrahlung. The external bremsstrahlung is also different, resulting 
in different momenta. Furthermore the external bremsstrahlung causes the tracks to have a 
non-zero impact parameter, which introduces a bias on the beam-constrained momentum. 

m 



The simulation should take into account all the differences between electrons and muons, 
when the momentum scale determined by muons is transferred to the electron momentum. 
However, mishandling any of these differences in the simulation may cause a difference 
between the electron momentum scale and the muon momentum scale, causing a discrepancy 
between the Z mass and E/p methods. In principle, the electron momentum scale can be 



checked using electron tracks. However, as shown in Table PCVI| , the uncertainties are too 
large to allow us to have concrete conclusions. 

This section describes the differences between electron tracks and muon tracks, how the 
simulation treats them, and the size of possible biases. 

Internal Bremsstrahlung Distribution: "Internal" photons are photons which are pro- 
duced at the vertex in a radiative W — *• ev^ event (or Z —>■ e^e^j event). For Monte 
Carlo events with no external photons, we find that the average E/p between 0.9 and 
1.1 is 1.0039. Part of this shift above 1, 0.0014, is from cut biases, and the internal 
bremsstrahlung shifts the peak by 0.0025. The distribution we are using would have to 
be wrong by ~100% for our fitted energy scale to come out shifted enough to account 
for the discrepancy between the energy scale from Mz and E/p. 

• The generator that is used for E/p simulation in these studies (PHOTOS in 



§§^Note that no material effects are included for the muons from the W and Z decays because 
they are negligible 
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two-photon mode) has been compared to the calculation by Berends and Kleiss 
of Reference , and the two generators give similar energy-angle distributions. 

Laporta and Odorico argue that inclusion of multiple photon radiation from 
the final state electron may change the energy loss distribution of the electron 



relative to a single photon calculation, such as Berends and Kleiss. Reference |]5^ 
contains an algorithm to calculate the effect of a cascade of final state photons. By 
construction, this algorithm reduces to Berends and Kleiss for the case of single 
photon emission. Their algorithm is implemented for W decays. The Laporta 
and Odorico case has the mean E/p between 0.9 and 1.1 lower by 0.00033. This is 
not insignificant, but it is not nearly large enough to account for the discrepancy 
between the Mz and E/p methods. The statistical error on the Monte Carlo for 
this calculation is 0.00015. 



Baur, Keller, and Wackeroth have done a calculation of the W — > ez/7 pro- 
cess which includes radiation from the W propagator. We have received their 
calculation in the form of a Monte Carlo The Monte Carlo can implement 



their calculation, and it can also implement Berends and Kleiss. We run sepa- 
rately in each mode and implement some simple model of CEM clustering of the 



photons and measurement resolutions. We find that ||3J] produces a value for the 
mean of E/p between 0.9 and 1.1 that is 0.00023 lower than the Berends and 
Kleiss result. 

External Bremsstrahlung Distribution: The formula we are using for the photon en- 



ergy distribution was calculated in 1974 by Tsai |^0|]. This formula is still referenced 
in papers written today, but it is possible that the formula is unexpectedly breaking 
down at high energies. Evidence that it is not is given by the SLAC measurement of 
the Landau-Pomeranchuk-Migdal effect described below ||5^ . They measured the rate 
and energy distribution of bremsstrahlung of 25 GeV electrons incident on different 
targets. For all the targets, they measured some level of bremsstrahlung suppression 
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at low photon energies, as expected, but at higher photon energies, their measured 
distributions agreed well with the expectation from pO[] . 

Low Energy Bremsstrahlung Cutoff: Since the number of external photons diverges as 
1/ we only consider external photons above a certain energy. In particular, we only 
simulate photons above y = 0.1%, where y is the fraction of the electron energy taken 
up by the photon. However, we can integrate the total fraction of the electron energy 
that is carried by photons below the cutoff. The total fraction is y = 0.1% x 0.085, 
where 0.085 is an approximation of the effective number of radiation lengths seen by 
the electrons, including the CTC gas and wires. We expect this to affect the energy 
scale by less than 0.0001, which is a negligible amount. As a simple check we have 
increased the cutoff and we do not see any significant change in the fitted energy scale. 
A similar argument holds for the internal photons. 

Beam Constraint Biasing E/p: The beam constraint can bias tracks that have un- 
dergone external radiation (bremsstrahlung) before the CTC active volume. 
Bremsstrahlung causes the tracks to have a non-zero impact parameter which biases 
the beam-constrained momentum. The simulation follows the same procedure, and so 
we expect this bias to be reproduced. Two possibilities are considered. 

• The radial distribution of material may be wrong. 

The average radius of external radiation (including half the CTC gas) occurs at 
22.21 cm in the simulation. The bias depends on r^, and so the location of the 
material might be sensitive to the scale. As a check the simulation is run with 
all the material before the CTC gas placed in the beampipe, or with all placed 
in the CTC inner can. The material is scaled so that (Xo) is the same for both 
cases, ftaii for the beampipe case is higher than the CTC case by about 1% of 
itself. The average E/p from 0.9 to 1.1 is higher in the beampipe case than the 
CTC case by 0.0003. Both of these changes are small. Considering that these 
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are extreme cases for variations in the possible distributions of the material, the 
expected changes are negligible. 

• In the simulation, the correlation between curvature and impact pa- 
rameter mismeasurement may not be correct. 

This would cause the Monte Carlo to produce the wrong bias from the beam 
constraint. However, in the Monte Carlo, we use CTC wire hit patterns from the 
real W data to derive a covariance matrix to use in the beam constraint. We use 
the identical procedure that is used to beam constraint the real data. The results 
are insensitive to the cuts on Dq and to variations of the correlation. 

We also try setting the energy scale with the E/p distribution before the beam con- 
straint. We compare the Monte Carlo distribution to the data distribution. We get a 
result for the energy scale which is consistent with the beam constrained E/p result. 

Landau-Pomeranchuk-Migdal Effect. Multiple scattering of the electron can suppress 
the production of bremsstrahlung at low photon energies [Q. Qualitatively, if the elec- 
tron is disturbed while in the "formation zone" of the photon, the bremsstrahlung will 
be suppressed. The "formation zone" is appreciable for the low energy bremsstrahlung. 
(Similarly, the electron bending in a magnetic field can also suppress low energy pho- 
tons, but the CDF magnet is not strong enough for this to be significant.) SLAC has 
measured this effect for 25 GeV electrons. The suppression of bremsstrahlung depends 
on the density of the material and occurs below y ~ 0.01 for gold and y = 0.001 for 
carbon, where y is the fraction of the electron energy taken up by the photon. The 
average density of material in the CDF detector before the CTC is closer to carbon 
than gold, and since we have a cutoff at y = 0.001, we are in effect simulating 100% 
suppression for the carbon case. This is a negligible effect on E/p. Any effect, if there 
were, will make the discrepancy bigger. 

Synchrotron Radiation. We considered the possibility that secondary particles, such as 
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synchrotron photons, may interact in the drift chamber, generating spurious hits and 
biasing the electron momentum measurement. To estimate the effect of synchrotron 
photons, we used a simple Monte Carlo simulation to convolute the synchrotron ra- 
diation spectrum for 35 GeV electrons with the photoelectric absorption length in 
argon/ethane. Assuming each absorbed photon to produce one drift chamber hit 
(except for the merging of nearby hits due to finite pulse widths), electron and photo- 
electron hits were fed to a hit-level drift chamber simulation and processed by the full 
track reconstruction software. The predicted bias in beam-constrained momenta due 
to synchrotron photons was ~ —0.02%, more than an order of magnitude too small to 
explain the energy scale discrepancy. We performed a second study, using a GEANT- 
based detector simulation under development for a future run of the CDF experiment. 
We used GEANT to simulate secondary particles near a 35 GeV electron, using the 
material distribution of the upgraded detector, and transplanted the secondaries into 
the same hit-level simulation used in the first study. The bias due to secondary parti- 
cles was again ~ —0.02%. We conclude that interactions of secondary particles in the 
drift chamber are unlikely to be the source of the discrepancy. 

Significant Energy Loss in Silicon Crystals. An electron moving through the material 
before the CTC will pass through ^ 400 fim of aligned silicon crystals. If it travels 
through the crystal along a major axis of symmetry, it can potentially lose significantly 



more energy than is lost through bremsstrahlung |53|. However, in the data we do not 
see any significant difference between electrons that pass through the SVX' and those 
that do not, relative to the Monte Carlo. This indicates that this is not a significant 
effect. 

Track Quality Comparison. In a completely data-driven study, we examined a large 
number of track quality variables, such as hit residuals signed in various ways, track 
X^, and correlations between hit residuals, as well as occupancies and pulse widths. 
While we had no quantitative model in mind to set the scale for comparisons, none 
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of the track variables we considered showed any significant difference between the W 
electron and W muon samples. 

4. Other Checks 

Inveiriant Mass Measurement: Calculating the invariant mass of Z — > e'^e" events 
makes use of a different set of track parameters than calculating E/p, and one could 
hypothesize errors in the angular variables causing errors in the invariant mass. We 
would not necessarily expect the electron and muon invariant masses to look the same 
since one uses Et and the other px- One could also imagine measurement correlations 
between the different tracking parameters which have the net effect of shifting the 
measured mass. The two tracks themselves could also be correlated since for Z events 
they are largely back-to-back. For example, if one track enters a superlayer on the 
right side of a cell, the other track will be biased to do the same. However, we have 
not been able to see any effect on the Z mass in the data. 

Inner Superlayers: Wires of the CTC inner super layers have larger occupancy than those 
of the outer superlayers, giving a higher probability of using wrong hits in the inner 
layers. To check this the Z electron tracks are refit with superlayers and 1 removed. 
While the resolution becomes worse, no significant change is seen in the means of E/p 
of the electrons or the invariant mass of Z electron tracks. Refitting is also done with 
the same tracks but by removing superlayer 5 instead of and 1. Again no significant 
change was observed in the means of E/p, or the invariant mass of Z electron tracks. 
The mean of the E/p distribution of W data is checked with the number of stereo or 
axial hits used in the track reconstruction. It is found to be insensitive to the number 
of hits. 

Coding Errors. Several independent E/p simulation codes produce highly consistent re- 
sults. 
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CEM Non-Linearity. When we applied the non-hnearity correction of Section [IV E| , the 
CEM energy scale factor as determined from E/p moved from 0.9963 to 0.9948, which 
makes the discrepancy between E/p and Mz worse. The uncertainty on the energy 
scale was also significantly increased by the uncertainty on the non-linearity. If we 
do not consider a non-linearity correction, then the discrepancy between the Z mass 
energy scale and the E/p energy scale is closer to 3.3 standard deviations. The data 
(see Figure ^), however, support a CEM non-linearity. 

Amount of Material is Incorrect. To increase the fitted energy scale by 0.5 %, we would 
have to increase the amount of material in the Monte Carlo by ~ 5.6 % of a radiation 
length. However, the tail of the E/p distribution of the W data is not consistent 
with such an increase. Moreover, the low tail of the invariant mass distribution of 
J/ip e^e~ decays (see Figure |5^) has been examined, and such an increase in the 
amount of material would significantly contradict the data. 

Backgrounds are Biasing the Result. It is possible that our estimate of the E/p shape 
of the background is flawed, and that there is a significant source of non-electron 
background in the E/p peak region that is biasing our energy scale fit. We consider 
the worst case possibility that all the background is located at one of the edges of 
the E/p fit region. To increase the Se{E/p) to 1, we would need to have about 
6% background piled up at E/p = 1.1. This is a factor of ~17 larger than the QCD 
background we have measured, and since we expect the QCD background to be largely 
fiat in E/p, we do not expect that backgrounds are significantly biasing our result. 
The agreement of the Z E/p fit with the W fit also indicates that the backgrounds 
are not a significant effect in the W fit. 

Tracking Resolutions Not Simulated Correctly. For the Monte Carlo, we smear the 
track parameters according to the calculated covariance matrix, and we then apply 
the beam constraint according to this same covariance matrix. Thus, in the Monte 



132 



Carlo, the covariance matrix used in the beam constraint describes the correlations and 
resolutions of the track parameters exactly. On the other hand, it is not necessarily 
the case for the data that the correlations and resolutions are described correctly by 
the covariance matrix. 

We can measure the correlation between impact parameter and curvature by plotting 
the average of qDq as a function of E /p. The slope of this plot for the data is slightly 
different than for the Monte Carlo. Since the Monte Carlo covariance matrix is the 
same matrix that is used to beam constrain the data, we conclude that the beam 
constraint covariance matrix does not perfectly describe the underlying measurement 
correlations of the data. 

To see how much of an effect this has on E/p we run the Monte Carlo as follows: 
We smear the Monte Carlo according to an adjusted covariance matrix, where all the 
off-diagonal terms are set to except for cr^(C, Dq), and which we fix according to the 
W data. When we apply the beam constraint, however, we use the same covariance 
matrices that are used by the data to do the beam constraint. In this way, we simulate 
the data more closely: smearing according to one matrix, and beam constraining 
according to a slightly different matrix. We find no effect on the average E/p between 
0.9 and 1.1. 

The Solenoid May Cause Non-Linearity in Photon Response. The solenoid coil 
presents ~ 1 radiation length for electrons in W and Z events, and also for any 
associated soft photons. Electron energy losses in the solenoid are not expected to 
affect our results since they are part of the CEM scale, which we are fitting for. 
However, it is possible that the soft photons are not making it through the solenoid 
and that this is distorting the E/p shape. As a simple check, we use a formula from 
the PDG Full Listings which describes the energy loss profile of a particle as a 
function of its depth in radiation lengths. We apply this formula to all the photons 
created in the Monte Carlo and reduce their energy accordingly. This is not a rigorous 
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check since we are applying the formula to low energy photons, which are in an energy 
region where the formula is not necessarily accurate. We rerun the Z Monte Carlo 
with this effect put in, and we treat this new Monte Carlo as "data" and fit it with 
the default Monte Carlo. Fitting gives a Monte Carlo energy scale of 0.99960, 
and fitting Mz gives a scale of 0.99935. We are interested in Mz relative to -B/p, and 
thus 0.99960 - 0.99935 = 0.00025 ± 0.00015. This is more than an order of magnitude 
too small to explain the energy scale discrepancy. 

5. Conclusion 

We have measured the energy scale using the peak of the Ej^ distribution of W data. 
The Ej'p distribution of Z events gives consistent results for the distribution of W 
events. However, if we set the energy scale with E jp^ then the invariant mass distribution 
of the Z events comes out significantly low. As a check we have refit the Run lA data with 
the Run IB Monte Carlo simulation, and the result agrees excellently with the pubhshed 
results. 

We have discussed several possible reasons that the Z mass comes out wrong. The 
problem could be a momentum scale problem or otherwise a tracking problem; it could 
be related to our simulation of Ej-p as presented in this paper; or it could be something 
theoretically unexpected. None of the plausible explanations considered here appears to be 
capable of creating a discrepancy of the magnitude observed in the Run IB data sample, 
and the source for the inconsistency remains an open question. 

For the final W mass measurement reported in this paper, we have used the invariant 
mass of the Z — > e^e~ and Z — > events. In this way, we have separated our energy 

scale measurement from almost all questions associated with the method. 
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